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ABSTRACT 


This  study  investigates  isothermal  and  non-i sothermal  miscible 
displacement  in  a  porous  medium.  The  effects  of  several  parameters  on 
dispersion  are  studied,  including  the  effect  of  a  temperature  front  on 
the  process.  The  work  done  includes  a  series  of  experimental 

measurements  of  thermal  and  concentration  profiles  for  two  component 
displacement  and  numerical  simulations  of  several  experiments.  The 
development  of  the  numerical  simulation  program  is  described.  This 
simulation  is  unique  in  that  the  equations  are  linked  by  temperature.  A 
dimensional  scaling  study  of  the  one  dimensional  equations  and  of  the 
displacement  parameters  is  also  included. 
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INTRODUCTION 


Miscible  floods  have  been  used  quite  successfully  in  some  types  of 
conventional  oil  reservoirs  to  improve  the  recovery  of  oil.  In  these 
floods,  solvent  (high  pressure  gas,  propane  or  liquid  carbon  dioxide)  is 
used  to  di spl ace  oil . 

For  heavy  oil  recovery  the  use  of  a  solvent  driven  by  steam  or  by 
a  hot  water  bank  has  been  proposed.  This  method  has  been  studied  by  a 
number  of  investigators,  including  Farouq  Ali(l).  It  has  the 

characteristics  of  both  miscible  displacement  and  steam  flooding. 
Attempts  to  study  this  process  by  scale  models  must  scale  both  heat 
transfer  and  the  displacement  process.  This  leads  to  difficulties  in 
satisfying  the  relevant  criteria. 

The  problem  was  therefore  defined  as  follows: 

"A  porous  medium  is  saturated  with  a  hydrocarbon  fluid  which  is 
displaced  by  a  second  hydrocarbon  fluid.  The  displacement  front  moves 
at  a  constant  velocity.  The  displacement  front  is  accompanied  by  a 
thermal  front  moving  at  a  different  velocity.  The  pressure  drop  within 
the  porous  medium  is  a  function  of  the  viscosity  of  the  moving  fluid. 
This  viscosity  is  a  function  of  fluid  temperature  and  composition". 

The  problem  was  approached  by  two  methods: 

1.  An  experimental  approach,  in  which  thirty-six  experiments  in  three 
sand  packs  were  used  to  explore  the  effects  of  fluid  velocity, 
inlet  mobility  ratio,  inlet  temperature  gradient  and  sand  particle 
size  on  dispersion.  Of  particular  interest  was  the  interaction 
between  the  displacement  front  and  the  thermal  front. 

2.  A  numerical  simulation,  in  which  the  equations  for  diffusion,  heat 
transfer,  and  pressure  drop  were  solved  by  a  finite  difference 
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method. 

The  goal  of  the  numerical  model  was  to  simulate  the  interactions 
of  the  thermal  front  and  the  dispersion  front,  as  determined 
experimentally.  The  study  included  a  scaling  analysis  of  the  problem 
which  demonstrated  the  incompatibility  of  the  scaling  groups. 


I.  SURVEY  OF  THE  LITERATURE 


The  use  of  laboratory  models  to  study  displacement  processes  in  oil 
reservoirs  requires  that  the  models  be  properly  scaled.  The  scaling 

process  determines  the  dimensionless  groups  relevant  to  the  process 
under  study  and  adjusts  the  variables  of  the  model  so  that  the  values  of 
the  dimensionless  groups  for  the  field  and  the  model  are  equal. 

There  are  two  methods  for  deriving  the  relevant  dimensionless 
groups.  The  first  method  is  dimensional  analysis.  In  this  method, 
variables  are  repeatedly  combined  until  a  useful  set  of  dimensionless 

groups  is  obtained.  This  technique  only  requires  that  the  complete  set 

of  relevant  variables  are  known.  The  procedure  has  been  described  by 
Greenkorn  (2)  and  by  Geertsma,  Croes  and  Schwartz  (3).  The  second 
method  is  inspectional  analysis.  In  this  procedure,  the  equations 
governing  the  system  under  study  are  written  in  dimensionless  form  and 
the  resulting  dimensionless  parameters  are  combined  to  yield  a  set  of 
dimensionless  groups.  Geertsma,  Croes  and  Schwartz  (3)  have  described 
this  procedure  and  have  derived  scaling  groups  for  waterflooding,  hot 
waterflooding,  and  miscible  displacement.  The  results  of  dimensional 
and  inspectional  analyses  on  the  non-i sothermal  miscible  displacement 
problem  can  be  found  in  Appendix  A. 

The  results  of  these  analyses  show  that  the  problem  cannot  be 
scaled  satisfactorily  due  to  conflicting  scaling  criteria.  The  fluid 
flow  and  heat  transfer  groups  demand  that  the  model  flow  velocity  and 
permeability  be  higher  than  the  prototype  values,  but  the  dispersion 
scaling  group  demands  that  the  term  vdp  (flow  rate  x  particle 
diameter)  be  equal  to  the  prototype  value.  Since  the  miscible 

displacement-heat  transfer  problem  cannot  be  adequately  modelled  in  the 
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laboratory,  it  is  necessary  to  predict  field  performance  by  a  numerical 
model . 

The  effect  of  different  variables  (velocity,  permeabil i ty ,  temper¬ 
ature,  viscosity)  on  the  dispersion  coefficient  k<j  is  of  interest  to 
anyone  who  would  like  to  predict  the  performance  of  a  laboratory  or 
field  experiment.  Miscible  displacement  experiments  have  been  described 
by  a  number  of  investigators.  Brigham,  Reed  and  Dew(4)  reported  experi¬ 
ments  in  which  the  parameters  velocity,  length,  particle  size,  viscosity 
ratio  and  core  diameter  were  studied.  They  determined  that: 

1.  At  favourable  viscosity  ratios  ( M< 1 )  the  displacement  was 
adequately  described  by  a  theoretical  solution. 

2.  The  dispersion  coefficient  k^  changed  by  a  factor  of  5.7  when  the 
mobility  ratio  was  changed  by  the  same  factor  (from  0.175  to 
0.998).  This  implied  a  linear  relationship  between  mobility  ratio 
and  dispersion. 

3.  When  the  viscosity  ratio  became  unfavourable,  the  theoretical  error 
function  curve  could  no  longer  describe  the  recovery  profile. 

4.  Viscous  fingering  or  bypassing  of  the  displaced  oil  by  solvent  was 
evident  at  unfavourable  mobility  ratios. 

5.  At  low  rates  of  flow,  the  dispersion  coefficient  k^  was  found  to 
be  independent  of  velocity. 

6.  At  higher  flow  rates  the  dispersion  coefficient  was  a  function  of 
velocity  and  of  particle  diameter. 

These  results  form  the  basis  of  the  equation  k^  =  D0  +  ovdp. 

Pozzi  and  Blackwell (6)  derived  a  set  of  scaling  criteria  for 
isothermal  miscible  displacement.  Their  experiments  were  performed 
primarily  in  the  unfavourable  mobility  ratio  region  and  were  used  to 
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validate  the  scaling  criteria  developed  by  inspectional  analysis.  The 
regions  in  which  these  scaling  criteria  are  valid  were  also  defined. 

Blackwell (5)  reported  an  experimental  study  in  which  a  dimension¬ 
less  dispersion  number  k^/Dg  was  plotted  against  a  dimensionless 
velocity  group  vdp/D0.  From  this  plot  the  regions  in  which 
diffusion  dominates  (vdp/D0  <  0.4),  and  in  which  convective  mixing 
dominates  (vdp/D0  >  4)  were  determined. 

Van  der  Pol  (7)  reported  results  which  indicate  that  diffusion  is 
the  dominant  mechanism  under  field  conditions.  He  also  reported  that 
small  core  experiments  in  both  the  mixing-dominated  range  and  the 
diffusion-dominated  range  produced  results  which  agreed  well  with 
Blackwell's  work.  It  was  concluded  that  the  dispersion  coefficients 
determined  in  laboratory  core  models  would  make  scaling  of  dispersion 
difficult  or  impossible  in  the  laboratory. 

Experiments  on  dispersion  of  thermal  energy  were  reported  by  Green, 
Perry  and  Babcock(8).  They  presented  a  numerical  solution  to  the  heat 
transfer  problem.  This  solution  was  a  convection-diffusion  equation  and 
provided  valid  solutions  as  long  as  stable  flow  prevailed.  Stable  flow 
was  maintained  by  injecting  a  chilled  liquid  in  order  to  produce  a 
favourable  mobility  ratio. 

The  calculations  of  thermal  dispersion  by  Green  et  al (8)  show  that 
thermal  diffusion  may  occur  through  two  mechanisms.  At  low  rates, 
conduction  predominates  and  the  static  conductivity  of  the  system  is  the 
effective  diffusion  rate.  A  plot  of  vdp  vs  k^  indicated  that 
at  higher  rates  dispersion  begins  to  dominate  the  heat  transfer 


process.  This  transition  occurs  at  vdp/kh  -  0.005 


cm3  C 
J 


At 


- 


6 


rates  below  the  value  of  vdp/^  =  0*0025,  the  heat  transfer 
coefficient  k^  may  be  represented  by  a  constant.  At  higher  rates,  the 
conduction  coefficient  will  include  a  velocity  dependent  term. 

Perkins  and  Johnston(9)  reviewed  a  number  of  papers  on  dispersion 
and  miscible  displacement,  and  covered  the  topics  discussed  above.  They 
also  discussed  the  effect  of  an  immobile  gas  saturation  on  dispersion 
coefficients. 

In  order  to  test  the  validity  of  the  equations  and  coefficients 
used  to  describe  the  miscible  displacement  process,  it  is  necessary  to 
solve  these  equations  using  finite  difference  techniques.  A  number  of 
authors  (10,  11,  12,  13)  have  presented  techniques  for  dealing  with  the 
equations  describing  the  heat  and  mass  transfer.  Peaceman  and 
Rachford(lO)  used  a  combined  backward  and  central  difference  approxi¬ 
mation  and  solved  the  resulting  equations  by  an  iterative  technique 
using  an  implicit  difference  formulation.  Their  solution  also  included 
a  pressure  solution  for  two-dimensional  isothermal  incompressible  flow. 

Chaud'nari ( 11)  discussed  some  of  the  limitations  of  Peaceman  and 
Rachford's  method  and  suggested  a  scheme  which  would  eliminate  most  of 
the  errors  inherent  in  this  technique.  Chaudhari  used  an  explicit 
finite  difference  technique  in  which  the  numerical  error  was  cancelled. 

Laumbach(12)  presented  a  high-accuracy,  semi-implicit  technique 
known  as  the  Truncation  Cancellation  Procedure.  This  method  uses  a 
system  of  equations  at  three  spatial  levels  and  two  time  levels  which 
gives  rise  to  a  tridiagonal  matrix.  The  results  presented  for  this 
solution  system  indicate  that  good  results  may  be  obtained  even  when  a 
small  number  of  points  are  used.  A  further  advantage  of  this  method  is 
that  it  is  stable  for  large  At  values.  This  procedure  was  used  by 
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Awang(13)  to  simulate  one-dimensional  heat  and  mass  transfer  in  a  porous 
medi urn. 

Awang  performed  several  heated  injection  miscible  displacements  in 
order  to  examine  non-i sothermal  dispersion.  A  numerical  simulation  of 
the  convection-diffusion  equation  for  both  heat  and  mass  transfer  was 
performed.  The  results  of  these  experiments  showed  that  a  dispersion 
coefficient  as  a  function  of  vdp  only  could  not  be  used  to  model  a 
heated  miscible  displacement.  The  results  also  indicated  that  the  local 
viscosity  ratio  was  not  constant  throughout  the  experiment  and  was  a 
function  of  fluid  composition  and  temperature.  The  dispersion  obtained 
in  these  experiments  was  found  to  be  temperature  sensitive,  as  it 
varied  with  the  thermal  gradient  in  the  pack.  The  results  of  this  work 
combined  with  the  results  of  Brigham,  Reed  and  Dew(4)  suggest  that  the 
thermal  gradient  in  the  pack  may  produce  a  viscosity  ratio  variation 
which  influences  the  local  dispersion  coefficient  k<j* 

Leventhal ( 14)  described  a  fourth  degree  accurate  procedure  known  as 
the  Operator  Compact  Implicit  (OCI)  method.  This  method  achieves  high 
accuracy  without  resorting  to  the  large  number  of  points  required  to 
solve  each  equation  by  standard  fourth  degree  implicit  systems.  In  the 
OCI  system  the  operators  defining  the  function  at  the  i-1  node,  the  i 
node,  and  the  i+1  node  are  generated  using  only  these  three  points 
rather  than  the  five  points  normally  required  for  a  fourth  degree 
one-dimensional  system.  Solutions  to  the  Buckley-Leverett  displacement 
problem  were  presented,  which  demonstrated  a  solution  which  is  more 
accurate  than  the  solutions  obtained  by  a  conventional  finite  difference 
method  or  by  the  Galerkin  method. 

Mi 1 1 er ( 24)  has  discussed  extensively  the  effect  of  condensation 
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fronts  and  thermal  fronts  on  the  stability  of  a  moving  fluid  interface 
in  a  porous  medium.  His  work  on  the  stability  of  a  steam  front  involved 
the  derivation  of  a  stability  criterion  in  which  he  assumed  a 
displacement  front  which  contained  a  small  wave-like  perturbation.  The 
regions  in  which  the  perturbation  would  grow,  decay  or  remain  constant 
were  then  delineated  by  a  summation  of  the  forces  tending  to  cause 
instability  or  to  damp  the  perturbation.  Miller  concluded  that  a 
decrease  in  volume  of  fluids  at  a  condensation  front  or  a  thermal  front 
would  tend  to  favour  the  stability  of  the  front.  He  presented  the 
example  of  a  steam  front,  which  would  remain  stable  in  spite  of  the 
unfavourable  viscosity  ratio  between  water  and  steam.  His  analysis  also 
considered  the  effect  of  the  heat  released  by  condensation  on  the 
stability  of  a  thermal  front. 

Mil  1 er( 23)  also  examined  the  stability  of  moving  interfaces  due  to 
phase  transformation  or  mass  transfer.  The  effect  of  interfacial 
tension  on  interface  stability  was  examined  and  it  was  found  that 
interfacial  tension  had  a  stabilizing  effect  in  most  systems. 

Miller  and  Jain(22)  described  the  stability  criterion  for  another 
type  of  system  consisting  of  the  vaporization  or  condensation  of  a 
dilute  fog.  They  determined  that  the  gain  or  loss  of  volume  due  to  the 
thermal  energy  absorbed  or  released  by  the  process  was  the  factor  which 
stabilized  or  destabilized  the  process. 
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II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  apparatus,  as  illustrated  in  Figure  1,  was  similar  to  that 

employed  by  Awang(13)  and  by  Brigham  et  al.(4)  Modifications  specific 

to  this  project  are: 

1.  A  4:1  flow  multiplying  cylinder  is  used  in  the  system  in  order  to 
provide  flow  rates  beyond  those  available  from  a  Ruska  pump. 

2.  The  insulated  core  assembly  is  enclosed  in  a  pressure  vessel  in 

order  to  allow  for  the  use  of  a  confining  overburden  pressure 
during  tests  at  high  injection  pressures.  This  modification  allows 
for  the  operation  of  the  system  at  high  temperatures  (above  the  oil 
bubble  point)  or  at  low  temperatures  (where  high  fluid  viscosities 
produce  high  pressure  drops).  The  arrangement  had  the  advantage  of 
allowing  for  the  use  of  a  positive  net  overburden  on  the  core  at 

all  times.  The  overburden  could  be  increased  as  required  to  pre¬ 
vent  internal  pressure  from  expanding  the  sleeve  and  thereby 
changing  the  packing  within  the  core.  This  possibility  was  con¬ 
sidered  to  be  more  serious  than  the  changes  in  permeability  caused 
by  overburden  pressure  variations. 

3.  The  core  sleeve  is  fabricated  from  0.010"  (0.25  mm)  stainless  steel 

in  order  to  provide  a  sleeve  wall  with  a  low  thermal  conductivity 

and  to  produce  a  durable  core  assembly.  Appendix  D  contains  heat 
transfer  calculations  considered  in  the  design  of  this  sleeve 
system. 

4.  The  sleeve  is  provided  with  fittings  through  which  six 

thermocouples  are  inserted. 

5.  The  core  assembly  was  mounted  vertically  and  fluids  were  displaced 
downwards  so  that  gravity  would  favour  the  stability  of  the  fluid 
fronts. 
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Figure  1  Schematic  of  Experimental  Apparatus 


Clockwise  from  lower  right 


Ruska  pump 

4:1  displacing  cylinder 
Sample  collector 
Constant  temperature  bath 
Pressure  transducer 
Electric  heater  plug 
Data  logger 
Pressure  vessel 
Insul ation 
Thermocoupl e 
Core  assembly 
Ng  overburden  supply 
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The  core  assembly  is  7.93  cm  in  diameter  and  78.7  cm  long.  The 
sleeve  is  packed  with  80-120  ym  (20-30  mesh)  or  20-30  ym  (80-120  mesh) 
Ottawa  silica  sand  by  means  of  a  particle  distributor  as  described  by 
Wygal(20).  This  packing  technique  consists  in  slowly  adding  sand  to  a 
dry  core  by  means  of  a  particle  distributor  consisting  of  five  layers  of 
coarse  screens.  The  filled  core  was  sealed  with  an  end  cap  and  the  core 
was  ready  for  leak  testing  and  insulation. 

Prior  to  running  displacement  experiments,  the  porosity  of  the  core 
was  determined  by  weighing  the  core  dry,  evacuating  it,  and  saturating 
the  core  with  deaerated  distilled  water.  The  core  was  then  reweighed 
and  the  volumes  of  sand  and  water  were  calculated  assuming  densities  of 
2.645  gm/cm3  for  quartz  and  0.9983  gm/cm^  for  water  at  25°C.  The  final 
porosities  were  found  to  be  0.328  for  the  80-120  ym  (20-30  mesh)  sand 

pack,  0.368  for  pack  2  and  0.356  for  pack  3  (20-30  ym  sand  packs).  The 

differences  in  porosities  suggest  that  the  particle  distributor  is  more 
effective  with  large  size  particles.  This  was  confirmed  by  packing  some 
practice  cores.  The  cores  were  reweighed  after  drying  to  ensure  that 

all  water  was  removed  prior  to  saturating  with  oil.  The  cores  were 

insulated,  placed  in  the  pressure  vessel  and  evacuated  before  saturating 
wi th  oil . 

The  hydrocarbons  used  in  these  experiments  were  a  group  of 
commercially  available  oils,  the  Gulf  Harmony  series  of  hydraulic  oils. 
These  oils  ranged  in  viscosity  from  25.  mPa*s  to  256.  mPa*s  at  40°C. 
Varsol  was  used  where  a  viscosity  of  1.  mPa*s  was  desired.  A  complete 
summary  of  oil  properties  appears  in  Table  1. 

The  experiments  were  carried  out  as  follows:  The  core  was 

saturated  with  the  oil  to  be  used  as  the  displaced  fluid  by  injecting 
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Table  1 

Properties  of  Oils  Used 


Oil 

Ki nematic 

Vi  sco si ty 
at  40°C 
mmz/s( cs) 

Densi ty 
at  40°C 
gm/cm3 

Ki  nemati  c 

Vi scosi ty 
at  100°C 
mm2/s(cs) 

Densi  ty 
at  100°C 
gm/cm3 

Refracti ve 
Index 
at  25°C 

Harmony  44 

30.48 

0.847 

5.28 

0.803 

1.4696 

Harmony  53 

65.3 

0.857 

8.71 

0.818 

1.4759 

Harmony  87 

210.3 

0.868 

19.12 

0.832 

1.4816 

Harmony  111 

302.5 

0.875 

24.95 

0.841 

1.4891 

Varsol 

1. 

0.767 

- 

- 

1.4400 

Note:  Viscosity  values  are  manufacturer ' s  data.  Density  values 
are  measurements  made  on  a  Haake  densitometer. 
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this  oil  at  v  =  .02  cm/s  until  the  refractive  index  of  the  effluent 

matched  that  of  the  injected  oil.  Velocity  is  defined  as  flow 

rate/core  cross  section.  This  displacement  required  a  minimum  of  three 
pore  volumes  when  the  fluid  initially  in  place  was  more  viscous  than  the 
injected  fluid.  The  displacing  cylinder  was  then  loaded  with  the  second 
oil  used.  The  constant  temperature  bath  was  adjusted  to  the  desired 
injection  temperature  and  injection  was  started. 

Effluent  samples  were  taken  at  one  minute  intervals  for  favourable 
mobility  experiments  and  at  three  minute  intervals  for  unfavourable 
mobility  runs.  Sampling  was  done  at  proportionately  longer  or  shorter 
intervals  for  higher  or  lower  rate  displacements  to  maintain  a  sampling 
interval  of  one  percent  pore  volume  for  favourable  mobility  tests. 
Sampling  was  continued  until  1.25  pore  volumes  were  displaced  in  the 
favourable  mobility  ratio  experiments  and  until  1.5  pore  volumes  for  the 
unfavourable  mobility  ratio  experiments.  Heated  injection  fluid  or 
chilled  injection  fluid  experiments  were  run  to  three  pore  volumes  in 

order  to  log  temperature  profiles  for  the  pack. 

The  composition  of  the  effluent  was  determined  by  refractive  index 
measurements.  Concentration  values  as  a  percentage  of  displacing  fluid 
were  obtained  from  a  calibration  curve  of  refractive  index  vs 

concentration  (Appendix  F). 

A  total  of  36  tests  were  performed.  Of  these,  17  runs  were  carried 
out  on  a  80-120  ym  sand  pack  and  19  runs  on  20-30  ym  sand  packs.  Four 
of  the  tests  in  the  80-120  ym  pack  were  repeat  experiments.  The  tests 
in  the  coarse  pack  included  experiments  at  four  favourable  and  two 
unfavourable  mobility  ratios.  Also  included  were  three  tests  with  a 
heated  displacing  fluid  and  one  test  with  chilled  displacing  fluid.  One 
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test  at  a  low  flow  rate  was  also  performed. 

The  20-30  pm  packs  were  used  for  experiments  at  six  favourable  and 
one  unfavourable  viscosity  ratio.  These  experiments  include  seven 
heated  injection  displacements  and  one  chilled  injection  displacement. 
One  low  rate  and  one  very  low  rate  experiment  were  included.  A  set  of 
three  displacements  in  a  water-wet  pack  was  also  performed. 


III.  DISCUSSION  OF  RESULTS 

A.  Experimental  Data 

The  concentration  of  displacing  fluid  in  the  effluent  from  all 
experiments  was  plotted  vs  time  on  linear  axes  to  yield  smooth  curves 
which  characterize  the  miscible  displacement  process  at  a  favourable 
mobility  ratio.  These  plots  closely  resemble  the  error  function 
solution  described  by  Brigham(21)  for  the  miscible  displacement 
problem.  Several  such  plots  are  shown.  Figures  2,  3  and  4  illustrate 
the  displacement  curves  obtained  from  an  isothermal  experiment,  a 
heated  injection  displacement  and  an  isothermal  displacement  at  a  lower 
velocity,  respectively.  The  corresponding  numerical  solutions  are  also 
shown.  Note  that  the  curves  are  plotted  using  a  cubic  spline  to  fit  the 
data  points. 

In  order  to  determine  the  dispersion  coefficients  for  each 
experiment,  the  concentration  data  was  plotted  as  concentration  vs  a 
pore  volume  number  EllR  on  probability  paper  (Appendix  F).  Volume 

/T^ 

displaced  at  a  given  time  is  V  and  Vp  is  the  fluid  contained  in  one 

pore  volume.  The  slope  of  the  resulting  straight  line  was  used  to 

calculate  the  dispersion  coefficient  kd  and  a  mixing  coefficient  for 

each  experiment.  (kd  =  -).  The  pore  volume  number 

3.625  (J) 

at  90%  concentration  is  V90  and  V10  is  the  pore  volume  number  at  10% 
concentration. 

The  displacements  carried  out  at  unfavourable  mobility  ratios  were 
also  plotted  on  linear  axes.  These  plots  did  not  produce  a  smooth 

curve,  and  exhibited  early  breakthrough  and  an  irregular  profile 
characteri stic  of  viscous  fingering,  (Figure  5).  When  the  data  from 
these  experiments  were  plotted  on  probability  paper,  the  resulting 
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Figure  2  NUtlERICRL  SOLUTION  FOR  R5 
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Figure  3  NUMERICAL  SOLUTION  FOR  R9 
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Figure  4  NUMERICAL  SOLUTION  FOR  R11 
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Figure  5  UNFRVOURflBLE  MOBILITY  DISPLRCEMENT,  R13 
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values  of  k^  and  a  were  found  to  be  very  large  and  should  be 
considered  to  be  the  averaged  results  for  viscous  fingering  and  not  as 
the  result  of  dispersion  at  a  stable  front. 

The  k^  values  or  their  dimensionless  counterpart  a  were  plotted 
(i)  vs  velocity  and  (ii)  vs  inlet  mobility  ratio  to  illustrate  the 
effects  of  velocity,  inlet  temperature,  sand  pack  size  and  inlet 
mobility  ratio  on  the  dispersion  process.  Table  2  is  a  complete  summary 
of  the  experimental  data.  Breakthrough  is  defined  as  the  earliest  time 
that  displacing  fluid  was  detected  in  the  effluent. 

Velocity  Effects 

The  effect  of  velocity  on  the  dispersion  coefficient  is  illustrated 
in  Figure  6.  The  displacement  experiments  at  three  velocities  are 
represented  in  this  plot  of  dispersion  coefficient  vs  velocity.  The 
dispersion  coefficient  appears  to  be  proportional  to  the  velocity  for 

the  20-30  and  80-120  ym  packs.  This  linear  relationship  between 
velocity  and  dispersion  was  discussed  by  Perkins  and  Johnston(9).  The 
plot  also  illustrated  an  effect  of  grain  size  on  mixing.  This  effect  is 
probably  due  to  a  decrease  in  packing  efficiency  in  the  20-30  ym  pack  as 
discussed  by  Brigham,  Reed  and  Dew(4)  who  noted  that  finer  sand  packs 
were  less  homogeneous. 

The  dispersion  plot  was  extrapolated  to  zero  velocity  to  obtain  the 
diffusion  coefficient  for  the  20-30  ym  sand  pack.  For  Harmony  87 

displacing  Harmony  44  in  a  20-30  ym  pack,  D0  =  0.00005  cm2/s.  This 

value  is  in  agreement  with  the  diffusion  coefficient  obtained  by  using 
the  diffusion  formula  found  in  (15),  p.574,  D0  =  A y"*^,  where  A  is 

the  reference  diffusivity  for  a  viscosity  of  one  cp.  This  diffusion  is 
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Table  2 

Experimental  Data 


Run 

Ti 

(°C) 

Inlet  M 

kd 

(cmz/s) 

pore  vol 
at  brkthru 

a 

V 

(cm/s) 

Comments 

1 

27.4 

0.089 

Pack  1 
'OVO'0 1'3' 

0.89 

0.84 

.0052 

(F) 

80-120  pm 

favourabl e 

2 

23.5 

11.2 

1.09 

167. 

.0215 

(U) 

unfavourabl  e 

3 

92.0 

0.089 

0.00136 

0.87 

0.86 

.0052 

F 

heated  run 

4 

25.6 

11.22 

0.267 

0.29 

167. 

.0052 

U 

5 

27.2 

0.110 

0.00165 

0.86 

1.02 

.0053 

F 

87  disp.  44 

6 

26.7 

0.708 

0.00188 

0.85 

1.30 

.0048 

F 

111  disp.  87 

7 

27.5 

708.0 

2.28 

0.15 

1380 

.0054 

U 

8 

28.2 

0.0014 

0.0072 

0.95 

0.46 

.0054 

F 

9 

94.0 

0.110 

0.0014 

0.91 

0.87 

.0054 

F 

heated 

10 

96.0 

0.708 

0.00226 

0.94 

1.38 

.0054 

F 

heated 

11 

26.1 

0.110 

0.00032 

0.92 

0.81 

.00127 

F 

low  rate 

12 

19.5 

0.708 

0.0015 

0.91 

0.93 

.00538 

F 

chilled  run 

13 

90.0 

11.2 

0.84 

0.34 

507. 

.00547 

U 

heated 

14 

26.7 

0.110 

0.00134 

0.91 

0.800 

.00550 

F 

repeat  of  5 

15 

28.0 

0.708 

0.00284 

0.84 

1.71 

.0054 

F 

repeat  of  6 

16 

27.8 

0.089 

0.00132 

0.90 

0.81 

.0054 

F 

repeat  of  1 

17 

27.2 

0.708 

0.00214 

0.91 

1.34 

.0052 

F 

repeat  of  6 

18 

26.8 

0.110 

Pack  2 
0.0O0454 

0.94 

1.27 

.0052 

F 

20-30  pm 

19 

25.3 

0.708 

0.00079 

0.90 

2.32 

.0049 

F 

20 

24.4 

11.22 

0.180 

0.41 

518. 

.0050 

U 

21 

91.5 

0.110 

0.000539 

0.95 

1.41 

.0055 

F 

heated 

22 

92.4 

0.708 

0.00139 

0.90 

3.66 

.0055 

F 

heated 

23 

89.0 

11.2 

0.249 

0.35 

652. 

.0055 

U 

heated 

24 

26.9 

0.110 

0.000408 

0.92 

1.19 

.0050 

F 

c  hi 1 1 ed 

25 

30.0 

0.110 

0.000183 

0.92 

2.35 

.0011 

F 

low  rate 

26 

29.3 

0.110 

0.000053 

0.94 

2.16 

.00035 

F  very  low  rate 

27 

29.2 

0.110 

Pack  3 
0.00092 

0.93 

2.14 

.0062 

F 

20-30  pm 
water  wet 

28 

27.5 

0.708 

0.00163 

0.92 

3.96 

.0061 

F 

water  wet 

29 

27.8 

11.22 

0.21 

0.38 

498. 

.0061 

U 

water  wet 

30 

92.0 

0.028 

0.00072 

0.92 

1.76 

.0057 

F 

heated 

31 

25.8 

0.116 

0.00055 

0.93 

1.42 

.0054 

F 

repeat  of  18 

32 

27.5 

0.708 

0.00118 

0.92 

3.18 

.0053 

F 

repeat  of  19 

33 

86.4 

0.089 

0.00054 

0.94 

1.40 

.0054 

F 

34 

92.5 

0.375 

0.00072 

0.93 

1.78 

.0056 

F 

53  disp.  44 

35 

89.2 

0.290 

0.00050 

0.94 

1.28 

.0053 

F 

87  disp.  53 

36 

91.5 

0.028 

0.00046 

0.91 

1.14 

.0057 

F 

heated 

. 
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Figure  6  Effect  of  Velocity  on  Dispersion 


23 


too  small  to  have  an  effect  on  mixing  at  experimental  flow  rates. 

Mixing  at  an  Unfavourable  Yiscosity  Ratio 

The  effect  of  an  unfavourable  mobility  ratio  on  the  mixing 

coefficient  is  illustrated  in  Figure  7.  The  heated  injection 

experiments  tests  (13,  23)  display  more  mixing  than  the  matching 

isothermal  tests  (4,  20).  Here  the  mixing  coefficient  a  is  defined  as  a 
dimensionless  number  (a=k(j/vdp). 

Perkins  and  Johnston(9)  have  defined  a  as  k^  =  Cavdp  where  a  is 

a  measure  of  the  inhomogeneity  of  the  pack.  In  this  work  a  is  defined 

as  a  dimensionless  mixing  coefficient.  All  of  the  unfavourable  ratio 
displacements  exhibited  a  high  degree  of  mixing  and  early  breakthrough 
of  the  displacing  fluid  (Figure  5).  These  characteri sties  suggest 
viscous  fingering.  The  unfavourable  mobility  ratio  experiments  using  a 
heated  displacing  fluid  exhibited  higher  mixing  coefficients  than  the 
matching  isothermal  experiments. 

Dispersion  at  M=1.0 

In  the  region  of  favourable  mobility  ratios,  a  set  of  experiments 
was  performed  in  both  sand  packs.  The  results  of  these  tests  are 
plotted  in  Figure  8  to  illustrate  the  trend  of  increasing  mixing  as 
mobility  ratio  increased.  The  lines  through  the  two  sets  of  points  may 
be  extrapolated  to  provide  estimates  of  the  mixing  coefficients  as  the 
mobility  ratio  approaches  unity.  This  extrapolation  yields  a  value  of 
a=2.5  for  the  mixing  coefficient  in  the  20-30  ym  pack  and  a=1.5  for  the 
80-120  ym  pack.  These  values  are  of  the  same  magnitude  as  the  value  of 
1.75  for  sand  packs  suggested  by  B1 ackwel  1  ( 5) .  The  increase  of  a  with 
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Figure  7  Mixing  Coefficients  at  Unfavourable  Mobility  Ratio 


25 


o 


O 

-H 

■H 

ru 

cr 

■H 


Figure  8  Isothermal  Displacement  Data 
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decreasing  particle  size  agrees  with  the  results  discussed  by  Perkins 
and  Johnston(9)  who  noted  that  pack  inhomogeneity  tends  to  increase  as 
particle  size  is  reduced. 

The  20-30  ym  pack  was  found  to  have  a  small  residual  gas  saturation 
which  may  have  increased  the  amount  of  dispersion.  The  presence  of 
residual  gas  was  indicated  by  a  decrease  in  the  pack  pore  volume  as 
calculated  from  a  material  balance  on  some  experiments.  This  volume 
change  was  of  the  order  of  four  percent.  This  gas  would  dissolve  in  the 
oil  during  displacements  at  high  pressure,  and  would  cause  the  effluent 
from  the  core  to  appear  cloudy.  Reducing  the  overburden  pressure  on  the 
core  appeared  to  reduce  or  eliminate  the  problem.  The  gas  was  assumed 
to  be  nitrogen  from  the  overburden  which  entered  through  a  fitting 

leak.  Perkins  and  Johnston(9)  noted  that  a  residual  gas  saturation  had 
little  effect  on  dispersion,  so  this  problem  was  not  considered  serious. 

Effect  of  Inlet  Temperature 

The  effect  of  heating  or  cooling  the  inlet  fluid  was  also 

investigated.  A  plot  of  mixing  coefficient  vs  mobility  ratio  (Figure  9) 

for  sand  pack  2  illustrates  an  increase  in  mixing  when  a  heated  fluid 

was  injected.  A  decrease  in  mixing  was  also  noted  when  a  chilled  fluid 

was  injected.  The  magnitude  of  this  effect  was  quite  small  in  spite  of 

temperature  gradients  sufficient  to  produce  an  unfavourable  mobility 
ratio  at  the  inlet  to  the  core  during  the  heated  injection  fluid 

displacements.  Similar  behaviour  was  noted  in  sandpack  1.  The  heated 
injection  tests  in  sand  pack  3  did  not  appear  to  have  an  increased 
dispersion,  possibly  because  the  removal  of  residual  water  was  not 

completely  successful. 
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Figure  9  Effect  of  Inlet  Temperature  on  Dispersion 
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Movement  of  Thermal  Front 

The  temperature  profiles  at  several  points  in  the  core  were  plotted 
for  some  of  the  high  temperature  displacements  and  for  one  chilled  fluid 
displacement,  (Figures  10,  11,  12).  Tj_  is  defined  as  the  inlet  fluid 
temperature  and  T0  as  the  initial  core  temperature.  The  temperature 
front  moves  through  the  core  more  slowly  than  the  concentration  front  as 
may  be  seen  from  Figure  3  and  Figure  10.  Because  the  two  fronts  will  be 
coincident  only  at  the  start  of  the  displacement,  the  temperature  front 
can  only  affect  the  displacement  by  initiating  instabilities  in  the 
concentration  front  at  this  time.  At  later  stages  of  the  experiment  the 
thermal  front  will  lag  behind  the  concentration  front  so  the  two  fronts 
will  be  independent.  The  thermal  front  in  the  unfavourable  mobility 
displacements  did  not  differ  significantly  from  that  in  the  stable 
displacements  as  may  be  seen  from  Figures  10  and  11.  The  jogs  in  the 
curves  were  thought  to  be  errors  in  the  data  logging.  All  curves  were 
plotted  by  a  spine  fit. 

A  stability  analysis  of  the  thermal  and  concentration  fronts  was 
carried  out  by  Miller's  stability  criterion.  The  results  appearing  in 
Appendix  E  confirm  the  experimental  evidence  that  the  thermal  front  was 
stable  in  the  range  of  experiments  conducted.  A  stability  analysis 
of  the  concentration  front  was  performed  using  a  stability  criterion  for 
a  miscible  front  moving  vertically  downwards.  This  criterion  may  be 
obtained  from  Miller's  criterion  or  from  Chuoke ( 26 ) .  Since  the  fluids 
used  were  completely  miscible,  there  is  no  capillary  term  in  the 
stability  number.  The  resulting  equation  implies  that  sand  pack 

diameter  is  not  a  factor  in  the  stability  criterion.  This  is  not 
true  for  the  case  where  capillary  forces  are  present.  The  only 
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stabilizing  force  in  the  experimental  system  was  gravity,  and  the 
gravitational  force  was  calculated  to  be  much  smaller  than  the  viscous 
forces  present  at  the  experimental  flow  rates.  This  was  confirmed 
experimentally  by  the  finding  of  instabilities  in  all  the  unfavourable 
mobility  displacements  performed. 

Pressure  Profiles 

Favourable  mobility  isothermal  displacements  exhibited  a  continuous 
increase  of  pressure  with  time  until  the  core  was  completely  saturated 
with  the  more  viscous  fluid.  At  this  point  the  pressure  became 
constant.  Heated  displacements  displayed  the  same  pressure  increase, 
followed  by  a  pressure  decline  as  movement  of  the  thermal  front  reduced 
fluid  viscosity.  During  chilled  fluid  injections,  the  pressure  drop 
increased  as  the  low  temperature  front  moved  through  the  sand  pack. 

Dispersion  in  a  Water-Wet  Sand  Pack 

Three  experiments  were  performed  in  a  sand  pack  initially  saturated 
with  distilled  water  which  was  displaced  from  the  pack  by  the  oil  to  be 
used  as  the  displaced  phase.  Displacement  was  continued  until  the  water 
phase  was  reduced  to  an  immobile  saturation.  The  oil  in  the  core  was 
then  displaced  by  the  oil  to  be  used  as  the  displacing  fluid.  Effluent 
samples  -were  taken  at  intervals  and  a  displacement  profile  was 
calculated  from  the  effluent  compositions  as  measured  by  refractive 
index.  The  dispersion  coefficients  were  calculated  from  the 
concentration  vs  time  curves. 

The  experiments  in  water-wet  sand  packs  produced  results  similar  to 
those  for  single  phase  experiments.  The  dispersion  was  found  to 
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increase  as  the  mobility  ratio  became  less  favourable.  At  an 
unfavourable  mobility  ratio,  early  breakthrough  and  viscous  fingering 
were  observed.  The  displacement  tests  at  favourable  mobility  ratios 
were  found  to  have  a  higher  dispersion  than  similar  single  phase  oil 
displacements.  The  unfavourable  mobility  displacement  was  found  to  have 
a  lower  mixing  coefficient  than  the  matching  displacement  with  no 
residual  water  saturation.  It  is  possible  that  water  mobilized  by  the 
high  injection  rates  had  an  effect  on  this  test.  The  results  of 
tests  in  water  packs  27-29  agree  with  the  results  of  Kasraie(25)  who 
noted  that  dispersion  in  sand  packs  decreased,  then  increased  as  the 
water  saturation  in  the  core  was  increased  to  the  residual  saturation 
point.  The  effect  of  residual  saturations  was  not  explored  further  in 
this  work,  but  it  was  noted  that  pack  3  heated  injection  experiments  did 
not  behave  as  did  the  experiments  in  packs  1  and  2.  The  dispersion  in 
tests  33-35  appeared  to  be  less  than  the  corresponding  dispersion  for 
isothermal  tests.  This  suggested  that  the  removal  of  residual  water 
from  the  pack  with  isopropyl  alcohol  was  not  completely  successful. 


IV.  NUMERICAL  ANALYSIS 


A.  Numerical  Procedure 

A  major  portion  of  this  study  involved  a  numerical  simulation  of 
the  convection-diffusion  equations  for  simultaneous  heat  and  mass 
transfer,  when  a  liquid  at  a  constant  temperature  displaces  a  second 
miscible  liquid,  initially  at  a  different  temperature. 

The  equations  describing  the  process  are: 

Heat  Transfer 

3 /  k  aiT\  _  3  (pvcpT)  _  3  ( C(  1  -4) )  prcpr  +  <j>Pcp]T) 

3X \  h  3X  /  3X  Ft 

Mass  Transfer 

3  / k  3  (pc)  \  _  3  (pvC) 

3x  \  d  3x  )  dx 

Pressure 

pk  3p  \  =  +  3p 
u  sx  y  at 

or.  Velocity 

3_(pv)  =  _  3p_ 

3  X  3 1 


Gravity  has  been  neglected  as  it  is 
small  relative  to  the  pressure  forces 
at  experimental  flow  rates  . 


The  numerical  solution  technique  selected  was  the  truncation 
cancellation  procedure  described  by  Laumbach( 12) .  This  technique 
produces  a  finite  difference  approximation  to  the  convection-diffusion 
equation  which  solves  for  the  n+lst  time  step  using  a  semi -impl ici t 
formulation.  The  equations  are  solved  for  a  one-dimensional  system,  so 
the  solutions  are  valid  for  the  stable  flow  regime  only.  The  system  of 
equations  and  the  finite  difference  formulations  are  found  in  Appendix 
B.  Awang(13)  performed  a  truncation  error  analysis  on  this  system  of 
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equations  and  demonstrated  that  they  are  second  order  accurate  for  all 

values  of  ^  and  that  they  are  fourth  order  accurate  for  aj=(— +— ). 

6  3 

u)  is  a  weighting  term  for  the  truncation  error  cancellation.  The 

truncation  error  terms  and  their  application  to  the  convection  diffusion 
equation  appear  in  Appendix  B. 

The  computer  program  utilized  solved  for  temperature  and  displacing 
fluid  concentration  by  a  semi-imp! ici t  procedure.  The  velocity  values 
at  each  node  may  be  solved  directly,  or  the  pressures  in  the  system  may 
be  used  to  calculate  velocity.  The  temperature  and  concentration  at 

three  points  at  the  nth  time  level  are  used  to  solve  for  the  temperature 
and  concentration  at  the  n+lst  level.  A  tri diagonal  matrix  was 
generated  using  rock  and  fluid  properties  at  the  nth  time  level  and  the 
matrix  is  solved  by  the  Thomas  algorithm.  The  values  of  temperature, 
pressure  and  concentration  obtained  were  used  to  update  the  fluid  and 
rock  properties  and  the  results  were  used  to  update  the  solution 
matrix.  The  solution  was  repeated  until  a  convergence  criterion  of  AC 
=  .0001,  AT  =  .01  °C,  AP  =  .001  bar,  was  achieved,  (Appendix  F). 

A  one-dimensional  set  of  equations  was  established  using  the  system 
properties  calculated  at  each  time  step.  Solution  of  the  resulting 
matrix  yielded  temperature  and  concentration  values  for  the  next  time 
step.  The  fluid  flow  equation  was  solved  by  two  methods. 

1.  The  velocity  at  each  node  was  calculated  by  a  direct  solution 

approach. 

2.  The  pressure  equation  was  solved  by  an  implicit  procedure  and 

velocity  was  calculated  from  the  pressure  values.  A  description  of 
the  solution  method  appears  in  Appendix  B.  The  three  equations 
were  coupled  by  temperature. 
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The  parameters  distance  step  ax  and  time  step  At  for  each  set  of 
computer  runs  were  determined  by  comparing  the  numerical  solution 
obtained  at  a  given  value  of  with  the  analytical  solution  at  the 
same  value  of  k^.  The  validation  of  the  numerical  solution  is 
described  more  fully  in  the  Discussion  of  Numerical  Results. 

B.  Numerical  Results 
Tests  of  the  Numerical  Program 

A  comparison  of  the  numerical  solution  with  the  analytical  solution 
is  presented  in  Figures  (13  -  16).  These  plots  show  that: 

1.  The  boundary  condition  C ( 0 , 0 ) =0  causes  the  numerically  calculated 
concentration  profile  to  lag  behind  the  profile  produced  by  the 
analytical  solution  (Figure  15).  An  overshoot  may  also  be  produced 
in  the  numerical  solution. 

2.  The  boundary  condition  C ( 0 , t ) =1  causes  the  numerically  calculated 
profile  to  lead  the  analytical  solution  (Figure  13). 

3.  The  use  of  a  smaller  distance  increment  causes  the  numerically  pro¬ 
duced  profile  to  follow  the  analytical  solution  more  closely. 
Conversely,  an  increase  in  the  distance  step  (Figure  14)  increases 
the  error  and  can  cause  an  overshoot  in  the  numerical  solution. 

4.  Changing  the  boundary  condition  to  C ( 0 , 0 ) =0 . 5 ,  C(0,t)=1.0  resulted 
in  a  numerically  generated  profile  which  matched  the  analytical 
solution  (Figure  16).  This  boundary  condition  eliminates  the  lag 
or  lead  of  the  numerical  solution  due  to  the  finite  width  of  the 
distance  step  ax  in  the  numerical  solution. 

One  simulation  was  attempted  using  the  Crank-Nichol son  method,  but 
this  system  proved  unsatisfactory  as  the  large  numerical  dispersion 
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overshadowed  the  actual  dispersion.  The  effect  of  numerical  dispersion 
on  solutions  to  the  convection-diffusion  equation  and  on  similar 
equations  is  discussed  by  Peaceman( 27) .  Peaceman  noted  that  the  con¬ 
vection  term  in  a  convection-diffusion  equation  will  produce  a  numerical 
dispersion  in  several  types  of  finite  difference  solutions.  This  dis¬ 
persion  may  appear  as  a  smearing  of  the  concentration  profile,  as  an 
overshoot  or  undershoot  behaviour,  or  as  instability  in  the  solution. 

Comparisons  of  simulations  carried  out  at  varying  values  of  At  and 
Ax  indicate  that  some  runs  were  unstable  and  that  there  was  an  overshoot 
in  other  runs.  A  good  fit  for  the  numerical  solution  occurs  between 
A t=180  s  and  ax=4  cm  (Figures  16)  and  At=120  s  (Figure  17).  The 
improvement  in  accuracy  obtained  by  iterative  updating  is  also  evident 
(Figure  18).  A  stability  analysis  of  the  truncation  cancellation  pro¬ 
cedure  by  the  Karplus  Criterion  (Appendix  C)  indicates  that  the  system 
will  be  stable  if  At<300  s,  if  Ax=4  cm  and  v=.0056  cm/s.  The  numerical 
runs  performed  indicate  that  the  actual  region  of  stability  is  At<240  s 
(see  Table  3).  This  stability  criterion  may  be  most  easily  met  by 
adjusting  the  time  step  to  an  appropriate  value  for  a  given  set  of  run 
parameters.  The  Karplus  criterion  is  not  completely  rigorous  as  it  does 
not  include  boundary  conditions.  This  omission  is  a  possible  reason  for 
the  discrepancy  in  stability  times. 

Results  of  Numerical  Simulations 

The  numerical  program  was  used  to  simulate  some  of  the  experimental 
displacements  (Figures  2-4,  Figure  19)  to  assess  the  ability  of  the 
numerical  system  described  to  predict  the  effects  of  mobility  ratio  and 
inlet  temperature  on  the  displacement  process.  The  results  of  these 
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Table  3 

Numerical  Simulation  of  Experiments 


Run  Mo. 

Inlet  M. 

a 

V 

cm/s 

Comments 

1 

0.089 

1.61 

0.00523 

3 

0.089 

1.60 

0.00521 

heated  disp. 

5 

0.110 

1.65 

0.00528 

6 

0.708 

2.01 

0.00475 

9 

0.110 

1.66 

0.00544 

heated  disp. 

8 

0.0014 

1.59 

0.00538 

numerical  overshoot 

10 

0.708 

2.01 

0.00538 

heated  disp. 

11  - 

0.110 

1.67 

0.00127 

12 

0.708 

1.82 

0.00538 

chilled  disp. 

8A 

0.0014 

1.56 

0.00538 

repeat  of  R8  at 

1 arger  At 

3A 

0.089 

1.30 

0.C0521 

pressure  solution 
for  R3 
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Figure  19  NUMERICAL  SOLUTION  FOR  R3 


46 


simulations  appear  in  Table  4. 

The  mixing  coefficients  obtained  from  the  numerical  curves  were 
plotted  against  mobility  ratio  for  a  set  of  simulations  of  the  80-120  ym 
pack  (Figure  20).  The  numerical  simulation  predicted  several  effects 
that  are  present  in  the  set  of  experiments  performed: 

1.  The  coupling  of  the  heat  transfer,  velocity,  and  dispersion 
equations  allowed  the  numerical  system  to  predict  the  changes  in 
velocity  due  to  temperature  changes. 

2.  Similar  results  may  be  obtained  by  employing  a  pressure  solution 
and  obtaining  velocities  from  the  pressure  values. 

3.  The  dispersion  equation  k^=vdpa  was  not  able  to  predict  changes 

in  dispersion  due  to  injection  of  a  heated  or  chilled  fluid  or  due 
to  changes  in  the  inlet  mobility  ratio.  It  was  necessary  to  supply 
a  mixing  coefficient  by  modifying  the  dispersion  equation  by 
including  a  mobility  ratio  term.  This  term  was  obtained  by 

multiplying  the  local  mobility  ratio  at  each  node  by  the  local 
dispersion  at  each  node.  The  local  dispersion  coefficient  becomes 
kd=vdpcrM.  A  description  of  the  local  mobility  ratio 
calculation  appears  in  Appendix  B. 

4.  The  use  of  a  linear  local  mobility  ratio  factor  in  the  dispersion 
equation  at  each  node  in  the  numerical  solution  was  capable  of 
reproducing  the  experimentally  observed  effects  for  isothermal 
simulations.  This  local  dispersion  coefficient  was  also  capable  of 
reproducing  the  effect  of  injecting  heated  or  chilled  fluids  in 
displacement  experiments. 

An  increase  in  mixing  was  predicted  when  a  heated  displacing  fluid 
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Tab! e  4 

Numerical  Stability  Tests 


Run  No. 

V 

(cm/s) 

At 

(s) 

AX 

(cm) 

Comments 

S12 

0.0025 

120 

8 

overshoot 

S13 

0.0056 

360 

4 

unstabl e 

S14 

0.0056 

300 

4 

unstabl e 

S15 

0.0056 

240 

4 

marginally  unstable 

S16 

0.0056 

220 

4 

stab! e 

S17 

0.0056 

180 

4 

stab! e 

S18 

0.0056 

90 

4 

stabl e 

R8 

0.0054 

180 

4 

stabl e 

R8a 

0.0054 

240 

4 

stabl  e 
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Figure  20  Numerical  Simulation  Data  for  20-30  Mesh  Pack 
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was  used  and  a  decrease  in  mixing  when  a  chilled  fluid  was  used. 
The  magnitude  of  this  effect  was  greater  when  the  inlet  viscosity 
ratio  approached  unity. 

6.  Isothermal  simulations  at  differing  velocities  predicted  mixing 
coefficients  to  be  large  relative  to  the  coefficients  expected  from 
diffusion  only,  as  the  range  of  velocity  studied  was  well  above  the 
velocity  values  where  diffusion  is  significant. 

7.  The  numerical  system  displayed  a  tendency  to  overshoot  when 
simulation  of  displacements  using  large  distance  steps  was  attempt¬ 
ed  (Figure  14).  Simulation  of  extremely  favourable  mobility 
displacements  tended  to  display  some  overshoot  behaviour  (Table  3, 
test  8).  This  was  remedied  by  modifying  the  time  step  used  (Figure 
21,  test  8A).  It  was  not  possible  to  completely  eliminate  the 
overshoot  from  all  simulations. 

8.  The  results  of  the  simulations,  when  plotted  on  probability  paper, 
produced  a  slightly  curved  line  rather  than  the  straight  line 
characteri Stic  of  the  error  function  solution  or  of  the  experimen¬ 
tal  data,  (Appendix  F).  The  curvature  made  it  difficult  to 
estimate  the  exact  slope  of  the  numerically  generated  concentration 
curve.  The  presence  of  this  curvature  was  assumed  to  be  the  result 
of  numerical  error.  This  curvature  was  found  to  be  much  smaller 
when  the  solution  method  was  changed  from  a  single  iteration  to 
multiple  iterations. 

9.  The  numerical  simulation  did  not  account  for  the  effect  of  packing 
or  porosity  on  the  dispersion  coefficient.  This  was  due  to  the  use 
of  a  fixed  value  of  1.75  for  the  mixing  coefficient  a.  Experiments 
on  a  number  of  different  packs  would  be  required  in  order  to 


" 

■ 


50 


CD 

O 

CO 


o 

o 

CO 

LO 


O 

O 

00 


o 

a 

o 

-'3- 


o 

o 

CM 

CO 


u 

aj 

cn 


o  QJ 
o  e 
■<r  -H 

CM  | — 


O 

O 

CO 


o 

o 

CO 


CD 


U0T^8J1U33U03 


Figure  2i  NUMERICAL  SOLUTION  FOR  R8 
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determine  the  relationship  between  porosity  and  dispersion. 

Perkins  and  Johnston(9)  discussed  the  results  of  several 
investigators  on  the  effect  of  packing.  It  was  reported  that  for 
finer  particles,  a  increased.  This  general  trend  was  found  in  the 
mixing  coefficients  calculated  for  the  two  pack  sizes  used. 

10.  The  numerical  simulation  could  not  reproduce  the  curve 

characteri stic  of  an  unfavourable  mobility  ratio  displacement 
(Figure  7).  The  use  of  a  local  M  value  to  modify  the  dispersion 
coefficient  produced  a  curve  similar  to  the  stable  displacement 
curve  (Figure  22).  To  model  unstable  displacement,  it  would  be 

necessary  to  use  a  two-  or  three-dimensional  model. 

11.  The  model  did  not  simulate  heat  loss  from  the  core,  as  a  heat  loss 

term  was  not  included  in  the  convection-diffusion  equation.  This 

system  of  equations  predicted  a  higher  rate  of  thermal  front 

advance  than  that  found  experimentally.  This  may  have  been  due  to 

the  lack  of  a  heat  loss  model  in  the  numerical  simulation. 

However,  the  model  was  capable  of  predicting  that  the  thermal  front 

(Figure  23)  would  advance  more  slowly  than  the  concentration  front 

(Figure  22).  The  numerical  simulation  could  also  be  used  to 

generate  profiles  of  the  local  viscosity  ratio  in  the  pack  at  any 

given  time  (Figure  24).  This  plot  was  constructed  by  calculating 

the  fluid  viscosity  as  a  function  of  composition  and  temperature  at 

each  node  on  the  numerical  grid.  The  local  viscosity  ratio  at  each 

node  was  then  calculated  as  it] _ £  .  The  profile  indicates  a  ratio 

v-j-h 

of  less  than  unity  at  the  stable  displacement  front  and  a  ratio  of 
greater  than  unity  at  the  thermal  front.  A  viscosity  ratio  profile 
for  an  isothermal  experiment  is  also  included  in  Figure  24. 
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Figure  22  ADVERSE  NOBILITY  RATIO  SIMULATION 
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Figure  23  NUMERICAL  TEMPERATURE  PROFILES,  R13 
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Figure  24  LOCAL  VISCOSITY  RATIO  PROFILE  AT  T=3600  S 
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12.  The  numerical  simulations  which  used  fluid  velocities  from  a 

pressure  converged  much  more  slowly  than  did  the  solutions  using  a 
direct  solution  for  velocity.  Pressure  solutions  required  eight  to 
eleven  iterations  to  converge  to  the  tolerances  in  C,  0.0001,  in  T, 
0.01,  and  p,  0.001  while  the  velocity  solutions  required  only  three 
iterations  to  achieve  the  same  tolerances.  This  difference  in 

speed  of  convergence  was  thought  to  be  due  to  the  ability  of  the 
direct  solution  to  provide  velocity  values  at  the  current  iteration 
of  the  solution.  The  calculation  of  velocity  values  from  pressure 
values  meant  that  the  velocity  values  available  were  those  from  the 
previous  iteration.  This  plus  the  increased  complexity  of  the 
solution  made  the  pressure  solution  slower  to  converge. 

13.  The  numerical  solution  for  pressure  produced  results  (Figure  25) 
similar  to  the  results  determined  experimental  ly  (Figure  26). 
The  experimental  pressure  profile  contains  a  step  at  the  pressure 
peak  not  found  in  the  numerical  solution.  This  step  is  due  to  the 
resistance  to  flow  of  the  end  cap  assembly.  This  small  increment 
of  flow  resistance  is  important  as  the  resistance  of  the  sand  pack 
is  very  small  when  a  80-120  yin  sand  is  used.  Figure  27  is  a 
comparison  of  the  concentration  profile  obtained  from  a  numerical 
solution  solving  for  pressure  with  an  experimental  concentration 
profile. 

C.  Results  of  Scaling  Study 

The  scaling  study  in  Appendix  A  indicates  a  need  to  scale  up  both  v 
and  k  as  the  model  is  scaled  down.  This  would  allow  for  modelling  of 
pressure,  gravity,  heat  transfer  and  viscous  forces.  However,  it  is 
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Figure  25  PRESSURE  SOLUTION  FOR  R3 
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Figure  26  Experimental  Pressure  Profile,  R3 
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Figure  27  NUflERICRL  SOLUTION  FOR  R3 
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necessary  to  scale  down  dp  as  v  is  increased  in  order  to  model 

dispersion.  However,  this  implies  that  k  must  be  scaled  down,  which 
conflicts  with  the  remaining  scaling  groups. 

An  example  of  the  application  of  scaling  groups  for  heat  transfer, 
viscous  forces  and  gravity  forces  is  as  follows:  A  1  ym^  (1  darcy),  78 

m  long  section  flowing  at  .152  m/d  is  modelled  by  a  100  ymz  sand  pack 

flowing  at  .0176  m/s.  The  prototype  would  be  in  a  velocity  region  where 
diffusion  dommi nates  (D0  ><jvdp)  while  the  model  would  be  dominated 

by  mixing  or  dispersion  (avdp  >D0).  The  sand  pack  would  be  able  to 
scale  gravity  forces,  pressure  forces  and  viscous  forces,  but  would  be 
unable  to  scale  dispersion.  Conversely,  a  model  correctly  scaled  for 
dispersion  would  not  model  gravity  or  viscous  forces.  The  experiments 
support  these  results.  The  very  low  rate  experiment  (R25)  produced  a 
low  dispersion  value.  Attempts  to  scale  pressure  forces  by  operating  at 
a  higher  velocity  in  a  larger  particle  pack  (R5)  produced  much  higher 
values  of  dispersion. 


' 
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V.  CONCLUSIONS 


A.  Experimental  and  Numerical  Results: 

1.  The  relationship  kc[=D0+avdp  as  described  by  Perkins  and 

Johnston(9)  was  found  to  be  valid  for  both  sand  packs  (Figure  6). 

2.  The  high  values  of  the  dispersion  coefficient  and  the  irregular 

displacement  profiles  occurring  in  the  experiments  at  M>1  (Figure 
5)  suggested  the  presence  of  viscous  fingering.  This  was  noted  by 
Brigham,  Reed  and  Dew(4)  and  other  studies. 

3.  The  successful  prediction  of  the  increase  in  dispersion  with  an 

increase  in  mobility  ratio  (Figure  20)  tends  to  support  the  use  of 
a  linear  local  viscosity  ratio  term  in  the  calculation  of  local 

dispersion  coefficients.  However,  this  does  not  imply  that  k<j 
for  the  entire  sand  pack  varies  linearly  with  mobility  ratio. 

4.  The  modification  of  dispersion  coefficients  by  a  mobility  ratio 

term  is  also  useful  for  predicting  the  influence  of  temperature 
gradients  on  dispersion  (Figure  20). 

5.  The  differing  rates  of  advance  of  the  thermal  front  and  the  dis¬ 
placement  front  were  predicted  by  the  numerical  model.  The 

numerical  model  did  not  contain  a  heat  loss  term  and  therefore 

could  not  predict  heat  losses  from  the  sand  pack.  Had  a  heat  loss 
term  been  included,  the  numerical  model  would  probably  have 

predicted  a  smaller  effect  of  injection  temperature  on  dispersion. 
The  lack  of  a  heat  loss  term  in  the  model  and  the  difficulty  in 

obtaining  adiabatic  conditions  in  the  experimental  displacements 

meant  that  it  was  not  possible  to  make  quantitative  measurements  of 
thermal  dispersion  in  the  sand  pack. 

6.  The  regions  of  stability  for  the  numerical  method  used  were 
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evaluated  (Table  3)  and  were  found  to  be  dependent  on  the  fluid 
velocity  used. 

7.  The  use  of  a  pressure  solution  to  calculate  velocity  produced  a 
workable  solution,  though  convergence  of  this  numerical  system  was 
much  slower  than  convergence  for  a  system  using  a  direct  velocity 
sol ution. 

8.  The  pressure  solution  is  more  versatile  as  it  may  be  written  for 
slightly  compressible  flow,  and  could  be  extended  to  a  compressible 
flow  case,  while  the  direct  velocity  solution  is  limited  to 
one-dimensional  incompressible  flow. 

9.  While  the  numerical  system  used  contained  several  simplifying 

approximations,  it  was  sufficiently  accurate  to  allow  qualitative 
predictions  of  rate,  temperature  and  mobility  ratio  effects. 

B.  Scaling  Study 

The  implications  of  the  velocity  dependence  of  kj  and  the  scaling 
analysis  are  as  follows: 

From  Appendix  A,  the  value  of  k^  modeled  field*  ^  *<d=Do 
+avdp,  then  dparp-jc-|e  must  be  decreased  as  vm0(je-|  increases.  This 

generates  a  conflict  as  the  remaining  scaling  groups  require  that  both  v 
and  dp  be  increased  when  scaling  down  from  field  size  to  laboratory- 
size.  Therefore,  a  model  of  a  process  in  which  dispersion  is 
significant  will  be  limited  by  the  dependence  of  dispersion  on  flow 
rate.  Such  a  system  may  be  scaled  only  if  it  is  not  necessary  to  model 
heat  transfer,  gravity  forces,  or  viscous  forces.  Conversely,  a  system 
modelling  heat  transfer,  viscous  and  gravity  forces  will  not  be  able  to 
simultaneously  model  dispersion. 
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C.  Practical  Application 

1.  The  criterion  modeled  field  W1'^  se1:  a  sca1in9  limit  of  2:1 

if  the  "field"  system  is  Harmony  87  displacing  Harmony  44  in  a 
20-30  pm  (80-120  mesh)  sand  pack  displaced  at  0.00035  cm/s  and  the 
maximum  permissible  error  is  15%.  For  a  fluid  pair  with  a  higher 
diffusion  coefficient,  the  scaling  constraint  will  be  somewhat  less 
severe.  At  Do=0.001  cm2/s,  the  system  may  be  scaled  5:1  to 
maintain  the  same  error  (15%).  This  scaling  factor  is  still  insuf¬ 
ficient  for  a  convenient  laboratory  model  of  a  petroleum  reservoir. 

2.  The  differing  rates  of  movement  of  the  concentration  and  thermal 
fronts  and  the  resulting  lack  of  a  strong  temperature  effect  on 
dispersion  suggest  that  the  direct  injection  of  a  heated  solvent 
will  leave  the  thermal  front  behind  the  solvent-oil  interface.  The 
resulting  process  will  be  equivalent  to  a  cold  solvent  injection. 

3.  The  numerical  solution  of  the  heat  transfer-di ffusion-pressure 
problem  is  of  interest  to  several  types  of  reservoir  problems.  The 
solution  of  the  convection-diffusion  and  pressure  equations  in  two 
or  three  dimensions  would  allow  for  simulation  of  systems  with 
unfavourable  mobility  ratios  or  simulation  of  heterogeneous 
reservoirs.  The  simulation  of  thermal  recovery  schemes  where  there 
is  significant  miscibility  of  components  would  require  the  solution 
of  all  three  equations  in  two  or  three  dimensions. 

D.  Recommendations  For  Future  Work 

1.  The  development  of  a  two-or  three-dimensional  model  of  the  non- 
isothermal  displacement  problem  would  be  an  asset  to  a  thermal 
recovery  reservoir  simulation. 
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2.  An  investigation  into  the  fundamentals  of  the  dispersion  process 
using  an  analytical  or  numerical  model  of  a  reservoir  pore  space 
would  better  define  the  relationship  between  mobility  ratio  and 
di spersion. 

3.  The  difficulties  in  reproducing  experiments  in  this  work  imply  that 
the  effects  of  packing,  porosity,  and  residual  saturation  of  a 
second  phase,  are  parameters  which  must  be  carefully  controlled  in 
a  displacement  experiment. 

4.  The  measurement  of  thermal  diffusivity  in  a  sand  pack  will  require 
one  of  two  appraches: 

i.  Run  experiments  in  a  low  pressure  model  with  vacuum 
insul ation. 

ii.  Match  the  thermal  profiles  of  experiments  using  a  numerical 
model  containing  a  heat  loss  term. 
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VII.  APPENDIX  A 


A.  Inspectional  Analysis  of  Miscible  Displacement 

i.  Differential  Equations 

1.  Fluid  Displacement  Equation 


3_  /  k  a  (pC)  \  _  3_(  vPC) 

3x  \  d  3x  j  3 x 

2.  Heat  Transfer  Equation 

3 _ /  k  3T  \  _  3 _ (pvcpT ) 

3x  y  h  3 x  J  3x 

3.  Darcy's  Law 

*  2>L  -  -  L  (12.  *  ipg\ 

3 1  U  l  3X  J 

4.  Dispersion  Equation 


k*-°vdp  +  w 


i  i . 

Boundary 

Condi tions : 

1. 

T(x,0) 

=T0 

for 

all 

X 

2. 

T(0,t) 

=Ti 

for 

all 

t  >  0 

3. 

T(x,») 

=Ti 

for 

all 

X 

4. 

O 

X 

o 

=0 

for 

all 

X 

5. 

C(0,t) 

=Ci 

for 

all 

t  >  0 

6. 

C  (  X,oo) 

<1 

for 

all 

X 

7. 

p(L,t) 

=0 

for 

all 

t 

iii.  Define  Dimensionless  Variables 

T*  =  T~To 
Ti-To 


_  3  ( 4>p  C ) 

Ft 


=  3 _  C((l-<j>)prcpr  +  4>pcp)T] 

3 1 


where  4>  _  is  the  fluid  velocity  v 
3 1 
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p 

P*  =  —  where  X0  and  t0  are  arbitrary  values 

Pi 

iv.  Substitute  Dimensionless  Variables  into  Equations  and  Simplify: 


3  1 

f  kd  3 (pC*)  ^ _ 

3  (vpC*) 

3  ( 4>pC* ) 

3eX0  ' 

y  3eX0 

) 

3eX0 

3xt0 

3C* 

=  ^o^d 

32C* 

 tov 

ac* 

3i 

CM 

O 

X 

-©- 

3e2 

4)X0 

3e 

3  1 

A.  3 

(T*(Tl- 

■ 

—i 

o 

+ 

H 

o 

V  3 

(pcpv(T*(T  i-T0)+T0) ) 

3eX0  1 

v  h  3eX 

0 

y 

)  3eXo 

_  3 

[T*(T-T 

0)  ( (l-(j>)prcpr 

+  4>pcp)  ] 

3t  t 

0 

3T* 

 kht0 

l 

32T* 

3t 

X  2 

Ao 

cd-4,: 

)  Pp  Cpr 

+  4>pcp] 

d£2 

-  V 

pep 

3T* 

xo 

(1-4,) 

p p  epp 

+  4>pcp 

3  £ 

3.  T(e  ,0)  =  T*A+B  =  T0 

4.  T(0,t)  =  T*A+B  =  Ti 

5.  C( e  ,0)  =  C*  =  0 

6.  C(0,t)  =  C*  =  1 

7.  T(e,«)  =  T*A+B  =  Ti 

8.  C(£,°°)  =  C*  =  1 

.  i/LP! Pi  -Apg\ 

9t  t0  U  i  J 

3e  _  k  t0Pi  3P*  _  kApgt0 

2 

3x  u  X0  4>  uxo  $ 


9. 


■ 

■ 


■ 


or 
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lEl  =  +  Apg  X0 

3e  k  Pi  pi 

10.  a  =  lid_  _  5p_ 

vdp  vdp 

v.  Equate  all  Dimensionless  Parameters  to  Unity: 


1. 


<j)X0 


1 


2.  hi  =  i 
4>x0 


3.  tokh  _ \ _ 

X02  ( 4>p cp  +  ( l-<f> )  pr  cpr) 

4.  tjjv  _ pep _ 

X0  ( (!-<(>)  pr  cpr  +  4>P cp ) 


5.  k  topi  =  1 
u  X0  (j> 


6 .  k  Apg  tQ  _  ^ 

y  X04> 


7. 


kd 

V  dr 


8. 


V  dr 


=  1 


vi.  Results: 

If  we  assume  that  X0=L 


t  -  <M-‘ 

0  *7 


1. 
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Tl 


= 

4>l2 


2.  t0  = 


iL 

V 


I2  = 


tv 


?  t  _  L2  ( ( 1  -4> )  pr  cpr  +  <j>Pcp) 


T3 


=  tkh 


'0 


T4  = 


L2 

((1-4) 

P  f  ^ 

■  4pcp) 

xo 

((1-4) 

Pj»  cp^  + 

4pcp) 

V 

pep 

tv 

pep 

L 

TTTTT 

Pr  cpr  + 

4pcp) 

5. 


t  =  nL2<|> 

0  kZT 


tkAp 
yL  2  4 


6.  t, 


 ul-4 


kApg 


tkApg 

t6  = 

yL(j> 


If  we  wish  to  solve  for  x: 
1.  combine  1  and  2: 


v  = 

Ao  — 


V 


combine 

1  and  4: 

X0  =  kd 

(1-<|>)  pr  cpr  + 

4>pcP 

4)  v 

pep 

e2  = 

<P  pep 

kd 

( l-4> )  pr  cpr  + 

4pcp 

combine 

1  and  6: 

X0  = 

kApg 

Ko  -  XkApg 
kdP 

combine  3  and  2: 

X  =  ^kh  _ ]_ _ 

v  (l-<j>)  pr  cpr  +  4>Pcp 

Xv  4>pcp  +  ( 1  -4> )  pr  cpr 

4 

4: 


6: 

_  _ 

Apgk  (l-<f>)  pr  cpr  +  4>p c p 

XApgk  (1-<j>)  pr  cpr  +  d>pcp 
k^y  <j> 


kh 


combine  3  and 


x0  =  h.  i_ 

V  pep 


Xv  pep 

kh  r- 


combine  3  and 


y  =  khy 


^6 


■ 


I 
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7.  combine  5  and  2: 


X 


o 


= 

U  V 


_  _  Xvy 

e7  “  - 

kAp 

8. 

combine  5 

and  4: 

Y  -  kAp 

4>p c p  (l_cf))  p |r>  cpp 

0  yv 

4>pcp 

e8  =  XV^ 

fpcp 

kAp 

4>p c p  +  pr  cpr 

9. 

combine  5 

and  6: 

X  =  AP 
o 

Apg 

p  _  ApgX 

£g  -  - 

Ap 

v  i  i . 

Scaling  Results  from  Groups  Obtained: 

Let  L  be 

Scaled  by  a  Factor  of  A 

Let  X=L 

1. 

from  Tg , 

JifL  =  A2  (assume  k^2  =  k^ ) 
ti 

2. 

from  t2, 

v2  _  1 

V],  A 

3. 

from  n , 

kd2  =  i 
kdi 

4. 

from  t6 , 

^1=1  (assume  M1  =  y2,  Api 
ki  A 

Ap2) 


■ 
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5. 

from 

Ti>, 

Ap2 

=  A 

Api 

6. 

from 

£9  , 

Ap2 

=  A 

(assume 

Api  =  Ap2 ) 

A  Pi 

7. 

from 

e6  » 

k2_ 

=  1 

(assume 

Api  =  Ap2 ) 

ki 

A 

8. 

from 

e7  , 

V2_ 

=  1 

(assume 

khi  =  kh2) 

Vl 

A 

9. 

from 

£3  , 

kd2 

=  1 

kdi 

10. 

from 

!<d 

dp2 

=  A 

vdn 

dPi 

11.  from  Do  ,  =  A 

vdp  dPi 

Note: 

1.  Group  9  conflicts  with  the  other  scaling  groups  as  illustrated  by 
groups  10  and  11. 

2.  The  equation  used  is  one-dimensional,  therefore  the  geometric 
scaling  groups  do  not  appear. 

B.  Dimensional  Analysis  of  Non- Isothermal  Miscible  Displacement 
i.  Variables 


L 

1 ength 

cm 

d 

diameter 

cm 

y 

vi sco si ty 

gm/cm  s 

cp 

fluid  specific  heat 

J/gm  K 

P 

fluid  density 

gm/cm3 

’ 
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kd 

dispersion  coefficient 

cm2/s 

V 

fluid  velocity 

cm/s 

Ap 

pressure  drop 

gm/cm  s2 

AT 

temperature  differential 

K 

g 

gravi ty 

cm/s2 

k 

permeabil ity 

2 

urn 

kh 

matrix  thermal  conductivity 

J/cm  K  s 

ii.  Matrix  of  Variables 


iii.  Determinant  of  Matrix  to  Far  Right  of  Matrix  of  Variables 
0  0  0  1 

110-3 
0-100 
0  0  10 
D=1 

iv.  Number  of  Dimensionless  Groups  Obtained 
N=n-r 

where 

n  =  number  of  independant  variables 


' 

' 


75 


r  =  number  of  dimensions  required  to  express  them 
N= 12 -4=8 

v.  Coefficient  Equations 
M  :  A2+A7+A8+A12  =  0 

L  :  A*  -A2+A3+2A4+2A5+2A6+A7-A8+A9+A1q-3A12  =  0 
t  :  -A2+2A3-2A4-A5-3A7-2A8-A10  =  0 
AT:  -A4-A7-A11  =  0 

vi.  Solution  of  Coefficient  Equations 
A12  =  -A2-A7-A8 

A11  =  A7+A4 

A10  =  -A2-2A3 -2A4 -A5-3A7-2A8 
Ag  =  -Ai-A2+A3_A5-2Ae-A7 

vii.  Matrix  of  Solutions 


1 

2 

3 

4  5 

6 

7  8  9 

10 

11 

12 

L 

y 

g 

cp  kd 

k 

kh  AP  d 

V 

at 

p 

1 

1 

-1 

0 

0 

0 

2 

1 

-1 

-1 

0 

-1 

3 

1 

1 

-2 

0 

0 

4 

1 

0 

-2 

1 

0 

5 

1 

-1 

-1 

0 

0 

6 

1 

-2 

0 

0 

0 

7 

1  -1 

-3 

1 

-1 

8 

1  0 

-2 

0 

-1 

vi ii .Dimension! ess  Groups  Obtained 


d 


' 
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2.  1L_ 

vdp 

3.  9d 

v2 

4  cpAT 

5.  M 
dv 

6.  iL 

d2 

7 .  khAT 
pdv3 

8.  A£_ 

pv2 

ix.  Modify  Groups  to  Obtain  Meaningful  Results 
Assume : 

P  2  _  CP2  =  92_  =  = 

Pi  cpi  gi  k^i 

Assume  that  kinetic  energy  is  negligible,  therefore  groups  must  be 
arranged  so  that  kinetic  energy  terms  do  not  appear  explicitly. 
Also,  since  we  desire  that  t^je]  be  smaller  than  t^-g-j^,  we  do 
not  wish  to  have  the  dimensionless  permeability  group  to  appear. 
Useful  groups  are: 

L 
d 

kd 
Lv 


2. 


' 
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3  ^h 

pcpvL 

4.  £P_ 

pgL 

5.  £i!L 
Vp 

6>  pcpAT 
Ap 

7. 

dv2 


Note:  If  heat  transfer  and  mass  transfer  are  not  strongly  coupled,  group 
6  may  be  ignored, 
x.  Results: 


let  k  =  A 


from  1  k  =  A 
di 


from  3  k.  =  k 
V]_  A 


from  2  kk  =  1 

kdi 


from  4  -_P_2_  =  1 

Api 


from  7^  =  1 
ki  A 


P2  = 
Pi 


from  5 


1 
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Fluid  groups  must  also  be  duplicated  if  two  fluids  are  used. 

usi  -  US2 
U01  U02 

P  si  =  PS2 
P01  PO  2 

CPsl  =  CPS2 
CPoi  CP0  2 


VIII.  APPENDIX  B 


Development  of  Finite  Difference  Equations  from  Differential  Equations 
A.  Convection-Diffusion  Equation 


3  /  ^  3(pC)\  _  3  (vpC)  _  4>3  ( pC ) 


3  X 


3  X 


3X  3 1 

may  be  modified  to  yield:  (ignore  cross  derivatives) 
2 


kd  £-3".C-  +  f  p3kd  +  2k^  3p(T) 


-  vP 


3  x 


3  X 


3x 


3C 

3X 


where 


<pp 


d£ 

3 1 


3  p ( T )  _  3p  3T 

3  X  3T  3X 

The  semi  implicit  finite  difference  formula  is: 

kdpjCT}  -  2Ct?+1+C[?^];+Cnl._1-  2C^C1+1) 
2ax2 


p3kd  +  2kd  3p.(Jl  -  vP 
3 x  3x  I  4ax 


i  /rn+ 1  rn+l,rn  rn  \ 

L_(ci+rci-i+ci+rci-i} 


.  *i(c!|+1-c!})  +e 

At 

The  truncation  error  from  a  Taylor  series  expansion  is 


e  = 


p£k  +  2k  3prn_  .  vp  \ 


y  3  x 

3x 

p3kd  + 

2kd  3p(T) 

3  x 

3  x 

-  Vp 


8  3x3 1" 

AX2  33C 


3x' 


-  p<(> 


At2  33C 


24  3t: 


rearrange : 
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■ 


80 


e  = 


ft P^Jid  +  2kri  3p^ 

U  3  X  d  3x 


-  vp 


VAt“ 
— 8  4> 


'p3kd  +  ,  3p(T) 

Jx~~  d  ^ 


vp 


4>aX' 

6v 


V2pAt2  \  3 3 C 


244 


3x23t 


This  is  similar  to  Laumbach's  truncation  error  term. 


Let 


e  =  BB 


■  3  3  C 


3X23 1 


or 


=  BB_(C^J-2C'J+1+C"+1-  C?.1+2C'J-C'j+1) 

AX2At 


let 


BB  =  BB' 


AX2At 


let 


BB  = 


-V2A t  _  (  p3'<d 


+  2k, 


244Ax' 


3  X 


3  P  ( T ) 

3X 


-  Vp 


VAt 


8<J>AX' 


-  ( p8kd 
l  3  x 


+  2kd  ?£.(II  -  vp  \  * 


3  X 


6vAt 


1  et 


AA  =  iL 


4>p 

At 


let 


CC  = 


.  ^p 


2ax" 
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1  et 


DD  = 


p3kd 


3  X 


2kd 


3p(T) 

3  X 


-  Vp 


4ax 


rewrite  the  equation: 

AA(Cn+1-Cn)  +  BB(Cn+1-2Cn+1+Cn+1-Cn  +2C 

i  i  i-1  i  i+1  i-1  i  i+1 


=  CC(Cn+1-2Cn+1+Cn+1+Cn  -2Cn+Cn  ) 
i-1  i  i+1  i-1  i  i+1 


+  DD(Cn+1-Cn+1+Cn  -Cn  ) 
i+1  i-1  i+1  i-1 


or 


ACn  +BCn+CCn  =  DCn+1+ECn+1+FCn+1 
i-1  i  i+1  i-1  i  i+1 

for  any  point  within  the  system. 

At  the  inlet  node,  the  system  is  supplied  by  a  source  of  strength 
vp ;  C-iniet  =  constant  =  Cx  where  C n -j e is  a  fluid  source.  If  C  is 
fixed  at  the  inlet  boundary,  the  equation  at  node  two  becomes: 

AA(Cn+1-Cn)+BB(-2Cn+1+Cn+1+2Cn-Cn) 

2  2  2  3  2  3 

=  CC(Cn+1-2Cn+1+Cn-2Cn+2C  )+DD(Cn+1-2C  +Cn) 

3  2  3  2  1  3  1  3 

Rewrite  the  equation: 

A'-C  +B 1 Cn+CCn  =  ECn+1+FCn+1 
1  2  3  2  3 

where 

A,  =  2DD-2CC  B,  =  -AA-2BB+2CC 


2CC-AA 


2CC-AA 
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At  the  outlet  boundary  the  system  is  closed  to  diffusion  but  is 
open  to  mass  transport  via  fluid  convection.  The  convection-diffusion 
equation  at  the  outlet  becomes: 

AA( Cn+1-Cn)+BB(2Cn+1-2Cn+1-2Cn  +2Cn) 

i  i  i  - 1  i  i  - 1  i 

=  CC(2Cn+1-2Cn+1+2Cn  -2Cn) 

i-1  i  i-1  i 

+  2DD(-Cn+1+Cn+1+Cn-Cn  ) 
i-1  i  i  i-1 

Now  let: 

A..  =  -2BB-2CC+2DD 
-AA+2CC 

R»  =  -AA+2BB+2CC-2DD 
-AA+2CC 

n..  =  -2BB+CC-2DD 
-AA+2CC 

p H  =  -AA+2BB-2CC+20D 

-AA+2CC 

Rewrite  the  equations  in  the  following  form: 

A"Cn  +B"Cn  =  D"Cn+1+E"Cn+1 
i-1  i  i-1  i 


Arrange  the  system  of  equations  in  matrix  form: 


A'  C  +  BCn  +  CCn 

1  2  3 

cn+1 

2 

E' 

F1  0  0 

A  Cn  +  BCn  +  CCn  = 

2  3  4 

cn+1 

3 

D 

E  F  0 

A  Cn  +  BCn  +  CCn  = 
i-1  i  i  + 1 

Cn+1 

i 

0 

D  E  F 

A"  Cn  +B"Cn  +0 
i-1  i 

Cn+i 

i 

0 

0  D"  E' 

A"  Cn 
i-1 


' 
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The  tri diagonal  matrix  thus  formed  may  be  solved  by  the  Thomas 
Algorithm.  A  similar  solution  applies  to  the  conduction-convection 
equation  for  the  heat  transfer  problem. 


B.  Heat  Transfer  Equation 

k  frA  _  3_(pCpvT)  _  9 _ ( (•  ( 1  -  <j> )  pr  cpr  +  <j»pcp)T) 

h  3x  )  3x  3 1 

may  be  modified  to  obtain: 

I,  92T  +  / 3_kh_  _pcPv\^X  =  3 _ ( ( 4>pcp+( l-<j>)  Pr  cpr)T) 

3x2  V  3 X  /  3X  3 1 

The  semi-impl ici t  finite  difference  approximation  is: 

J^h _ (T^tl-2T?+1+T?tl+Tq _!-2T9+T^+1 ) 

2ax2 

+  Hh  _  PCPV  1 _ (Tr?ti-T!?tl+T,?+i-T,?_1) 

3  x  4ax 

=  M(Tf: l-T?)  +  e 

At 


where 


[M]  =  ( 4>p cp  +  ( 1  -4> )  pr  cpr) 

The  truncation  error  term  obtained  from  a  Taylor  series  expansion 


e  = 


=  -  [M] 


At' 


33T 


axVt 


24  3 1; 


e  = 


+ 


-(vpcp)2  At' 
[M]  24~ 


3kh  _  vPcp\  At^  vpcp  \  3 3T 

/  8  [M] 


3  X 


3X23t 


At233T 
8  3X3 1' 


' 
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1  et 


e  =  BB1  111 


3x^3t" 


or 


no i  /-rii+l  9Tn+l ,xn+l  -T- n  ,9-rn  -T-n  \ 

=  _  (Ti-l“2Ti  +Ti+l"Ti-l+ZTi"Ti+lJ 

Ax2At 


1  et 


BB  = 


BB' 


AX2At 


let 


BB  = 


[M]AX2  +  At2 VpCp 


2.2 


3kh  -vpcp\  -(vpcp)*At 


vpcp  8  [M]  y\3x 


[M]  24 


1  et 


AA 


<fpcp  +  (l-cj>)  pr  cpr  _  [M] 


At 


At 


1  et 


let 


DO  = 

4a  x 

rewrite  the  equation  as  for  the  convection-diffusion  problem: 

ATn  +BTn+CTn  =  DTn+1+ETn+1+FTn+1 
i-1  i  i+1  i-1  i  i+1 


C. 


Fluid  Flow  Equation  in  One  Dimension 

The  mass  balance  for  the  total  system  is: 


■ 
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3  (pv)  =  _  4>3 p 
3  x  at 

This  may  be  rewritten  as  the  continuity  equation  for  one  dimensional 
non-i sothermal  flow: 

V3p  +  p3v  _  —  4>3  p 

3  X  3  X  3 1 


where  p  =  p(T) 

or 

3p(T)  = 

3p 

3T 

3  X 

3T 

37 

3p(T)  = 

3p 

3T 

3  t 

3T 

dt 

In  finite  difference  form: 


,  n+1  n+1.  .  n+1  n+1.  ,  n+1  n+1 

( V 1  +^+v i -h)  3p  (Tj-n-Tj.p  +  pn+l 

2  3T  2ax  i  ax 


,  n+1  n  \ 
_  43  P  (T-j  -T  -j ) 


3  T 


At 


This  may  be  solved  directly  if  velocity  is  given  at  node  1 


D.  Pressure  Equation 

Derivation  of  the  pressure  equation  for  non-i sothermal  slightly 
compressible  flow 
for  a  slightly  compressible  fluid: 


3  (pv)  _  -<ppC  3 p 
3  x  3 1 


where  C  is  the  fluid  compressibility 


for  a  non-i  sothermal  fluid  flow: 


. 

' 

■ 


86 


9  (p v)  =  -<j>3p 

3  X  3 1 


where  p  =  P(T) 
combine  the  equations: 

Mp(J)v)  .  /  3pm  +pC  ap) 
3x  Vat  at  / 

Darcy's  Law 


v  =  -k_  9JP 
y  3x 

Note:  The  gravity  term  is  neglected  as  gravity  forces  are  small 
compared  to  viscous  forces  at  experimental  flow  rates. 


—  /Ten  L  iPJ 

3X  ^  p  3  X/ 


The  change  in  density  due  to  pressure  is  very  small  at  experimental 
pressure  drops,  so  the  fluid  density  should  be  written  as  a  function  of 
temperature  only.  The  dependence  of  fluid  density  on  variables  other 
than  pressure  complicates  the  problem  in  that  it  requires  an  iterative 
solution,  as  an  analytical  solution  will  become  intractable. 

Define  transmi ssibil ity  as  follows: 


P  1  -h 


i-k 


(p  j-i+p j)  2k 

2  (ui-i+ui)  ax 


Place  pressure  equation  in  finite  difference  form: 


i+h 


Tn+h 


/  n+1  j+l\ 
(Pi+rPi  )  - 


p  1  - 


Tr 


i-b 


/nn+l  n+1  \ 
(Pi  “Pi -1 ) 


4>  A  X 


( (Pi 


n+1 


C  + 


9p 

df 


(T?+1 


At 


i 


' 
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Boundary  Conditions: 

The  boundary  conditions  common  to  both  heat  and  mass  transfer  are 
as  follows: 

At  the  inlet: 

C(0,t)  =  Cx  =  1  C(x,0)  =  0 

T(0,t)  =  Ti  T ( x , 0 )  =  T0 

v(0,t)  =  At  the  outlet: 

p  ( L ,  t)  =  0 

The  equations  used  for  the  fluid  and  rock  properties  are  as 
fol 1 ows: 

1.  Fluid  Density 

p  =  A+BT+CT2  HL 

cirr 

The  constants  A,B  and  C  were  determined  for  each  oil  used  by 
fitting  a  least  squares  curve  to  the  densities  measured  at  22°C, 
50°C  and  70°C  for  each  oil  used. 

2.  Density  of  Fluid  Mixtures 

_  1  gm 

m  Ci/pi  +  ( 1 -C i ) /p  2  cm3 

This  is  an  adaptation  of  a  specific  volume  formula  found  in  "The 
Properties  of  Gases  and  Liquids"  p.  86  (15). 

3.  Fluid  Viscosities 

p  =  A  exp  (1) 

T 

This  correlation  is  known  as  Andrade's  equation,  (Ref.  15)  p.  437. 
The  constants  A  and  B  were  determined  by  fitting  the  viscosities 
found  in  the  manufacturer' s  literature  to  the  equation. 


. 
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4.  Viscosities  of  Mixtures 

L  =  1000  in  20 
1  in  \i\-in  .0005 


L  m  ~  *-i  +  (1-Ci  )l-2 


um  =  .0005  exp 


1000  in  20 


L 


m 


The  correlation  used  by  Awang  (Cragoe's  equation)  and  a  logarithmic 
correlation  (Equation  19-8,  Ref.  15)  were  compared  with  viscosities 
measured  with  a  Nametre  viscometer  for  several  oil  mixtures. 
Cragoe's  correlation  (above)  was  considered  to  be  more  accurate, 
(Appendix  F). 

5.  Diffusivity  of  the  Bulk  Oil 


D0  =  AuQ 


cnf 


where  q  is  valued  between  -0.5  and  -1.0;  from  Ref.  15  p.  574. 

6.  Dispersion  Coefficient  in  Sand  Pack 


kd  =  Ro  +  1.75  vdD  — 
a  F4>  p  s 


This  correlation  was  proposed  by  B 1 ackwel  1(5).  Brigham  et  al  (4) 
noted  that  the  dispersion  coefficient  kd  appeared  to  be  a 
function  of  the  mobility  ratio.  Therefore  the  dispersion 
coefficient  was  defined  as  follows: 

kd  =  Bo  +  1.75  vd0M  — 
a  F4>  p  s 
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7.  Thermal  Conductivity  of  Sand  Pack 

kht  =  .0189  -  .000017T  _ - _ 

cm*  s*°C 

From  data  of  Anand,  et  al  (15).  The  properties  of  Boise  standstone 
were  assumed  to  be  similar  to  that  of  the  sand  pack.  The  value 
obtained  here  is  slightly  higher  than  the  value  obtained  for  tar 
sand  by  De s ton ( 17).  These  values  were  considered  useful  approx¬ 
imations  as  no  correlations  for  heat  transfer  in  unconsolidated 
sand  packs  saturated  with  refined  oils  were  available. 

8.  Thermal  Diffusivity  of  Sand  Pack 

kh  =  kht  +  .172S(vd p)1'4  — * 

K  cm*s*  C 


This  equation  was  proposed  by  Green  et  al  (8)  as  a  best  fit  of  the 
data  they  obtained  for  thermal  diffusion  in  an  unconsolidated  sand 
pack. 

9.  Mobility  Ratio 

This  is  defined  as  the  oil  viscosity  at  the  i+1/2  node  divided  by 
the  viscosity  at  the  i - 1/2  node.  This  ratio  can  be  approximated 
by: 

M  = 

Ui-l+Ui 

10.  Heat  Capacity  of  Oil 

The  "Chemical  Engineer's  Handbook" ( 18)  p .3-136  lists  the  following 
correlation  for  the  thermal  capacity  of  oils: 

cp  =  _LZ18  +  .00377  (T-15)  J— 

/ p is  gm  C 
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p i b  is  the  oil  density  at  a  reference  temperature,  (15°C). 

11.  Heat  Capacity  of  Rock 

The  specific  heat  of  quartz  can  be  determined  from  a  correlation 

found  in  Ref.  (19)  p.  f-68: 

cpr  =  .753  +  .0025T  — — 

gm°C 

12.  Rock  Density 

A  formula  for  the  density  of  quartz  may  be  obtained  by  combining 
data  from  Ref. (19)  p.  f-68  and  Ref. (18)  p.  3-20: 

Pr  =  2. 645-5 . 5X10~7T  90L 

cm3 
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IX.  APPENDIX  C 


A.  Stability  Analysis 

A  stability  analysis  of  the  Truncation  Cancellation  Procedure  via 
the  Karplus  Criterion  (28): 

1.  Write  the  equation  in  the  form  used  by  the  matrix  solution: 

AA(Cn+1-Cn)+BB(Cn+1-2Cn+1+Cn+1-Cn  +2Cn-Cn  ) 
i  i  i-1  i  i+1  i-1  i  i+1 

=  CC(Cn+1-2Cn+1+Cn+1+Cn  -2Cn+Cn  ) 
i-1  i  i+1  i-1  i  i+1 

+  DD(Cn+1-Cn+1+Cn  -Cn  ) 
i+1  i-1  i+1  i-1 

Rearrange 

AA(Cn+1-Cn)+BB(Cn+1-Cn)-2BB(Cn+1-Cn) 
i  i  i-1  i  i  i 

+BB(Cn+1-Cn)-BB(Cn  -Cn)-BB(Cn  -Cn) 
i+1  i  i-1  i  i+1  i 

-CC(Cn+1-Cn)+2CC(Cn+1-Cn)-CC(Cn+1-Cn) 
i-1  i  i  i  i+1  i 

-CC(Cn  -Cn)-CC(Cn  -Cn)-DD(Cn+1-Cn) 
i-1  i  i+1  i  i+1  i 

+DD( Cn+1-Cn ) -DD(Cn  -Cn)+DD(Cn  -Cn)  =  0 
i-1  i  i+1  i  i-1  i 

Terms  will  cancel  to  yield  the  stability  criterion: 

-AA+2BB+2CC  =  0 


When  the  isothermal  case  is  considered;  the  stability  criterion 
becomes: 

-2j)p  +  V2p&t  +  k(jp  -r  q 

3At  6(^AX2  AX2 


or 
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-2<+>p  +  v2pAt2  +  k^pAt  *  n 

2  2 
3  6<i>AX  AX 

For  a  typical  experiment: 

<j>  =  .328 

p  =  .87  gm/cm 

v  =  .00561  cm/s 

Ax  =  4.0  cm 

k,j  =  .0030  cm2/s 

-  .1902  +  .8696  x  10“6  At2  +  .000163  At  <  0.  The  system  is  stable 
for  At  <  300  s. 

Actual  instability  occurs  at  At  =  240  s. 


X.  APPENDIX  D 


A.  Heat  Transfer  Calculations  for  Miscible  Displacement  Test  Core 

Heat  Transfer  Through  Sand 

k^t  =  *865  W/m  °K(0.5  btu/hr  ft  °F)  (Green,  Perry  &  Babcock, 
Ref.  8) 

Qsand  =  Kht  A  AT) /d 
Qsand  =  *865  -00492  70/. 0792 

^sand  =  3.77  W  (7.13  btu/hr)  where  d  is  the  core  diameter  in  m 
Ksieeve  =  16.5  w/m  (9.5  btu/hr  ft  °F)  for  316  stainless  steel 
Qsieeve  =  16-5  *000113  70/. 0792 

^sleeve  =  1-64  W  (3.12  btu/hr)  for  a  .018"  wall  stainless  steel 
si eeve 

for  a  .010"  thick  stainless  sleeve,  Qsieeve  =  *91  w 
for  a  .010"  sleeve,  QS1 eeve^Qsand  =  -24 
The  .010"  sleeve  chosen  in  order  to  minimize  the  contribution  of 
the  sleeve  to  total  heat  conduction,  and  to  allow  the  overburden 
pressure  used  to  act  on  the  core. 
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Figure  28  Core  Heat  Transfer  Problem 


^sleeve 

Q 

xore 


must  be  small  for  experimental  accuracy 


XI.  APPENDIX  E 


A.  Miller  Stability  Criterion  for  Stable  Thermal  Displacement  Front 


JL  =  25°C 

Tz 

=  90°C 

Pi  =  .888  gm/cm3 

Pz 

=  .847  gm/cm3 

Ml  =  8.00  9™ 

M2 

=  .34  9™ 

cm  s 

cm  s 

Vx  =  .00540  cm/s 

V2 

=  .00561  cm/s 

ki  =  100  urn2  (darcys) 

k2 

=  100  urn2 

k^2  =  •°20  J/cm  s  °C 

v2 

=  .005  cm/s 

(relative  velocity  of  rock  to  thermal  front) 
cpi  =  2.01  J/gm  °C  cp2  =  2.25  J/grn  °C 

The  system  is  stable  when: 

(gravity  forces)  +  (viscous  forces)  -  (stabilizing  effect  of  heat 
transfer)  >  0.0 
Gravity  term 

-  ( p 2  -  pi)  g  a 
=  (.847  -  .888)  981.0  a 

=  +  40.18  aiE _ 

3  2 

cm  sec 

=  +  40.18  a  .^-.ll- 
Viscosity  term 

\  a  108 
k2  ! 

=  -41,200  a 

cirr 

Thermal  stability 
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yi 


U2 


P 1 k 1  p  2k2 


kh  tW2-y)  (Tl-T2) 

Hi-H2-cp2(Ti-T2)  (I2=X) 

a 


10 


8 


8.00  +  .34  \  .020  y(Y2-y)  (-65) 

.888x100  .847x100  )  -127.32-2.25  (-65)  [1211) 

'  a 


10* 


Y 


=  P2  cp2  /V-V2_\+  Pr  CP r  v 

kh  l  t  kh 


convective  heat  transfer 
conductive  heat  transfer 


Y  =  -.722 


Y2  =  I  +  A*  + 


Y 2  =  .5096  or  y 2  -  .175  for  x  -  2.  dsanc|pack  or  x  ^sandpack 

cl  =  JL  (X  is  the  largest  probable  wave  length) 

X 


a  =  .7923  or  a  =  .3961 
Gravity  forces 

-  ( P 2  -  Pi)  9  a 
=  +15.9  dynes/cm2 
Viscous  forces 


a  /  -  »11L  +  \  108 

ki  k2 


=  -  16310. 


dynes 


cm 


Stability  due  to  heat  transfer 


yi  +  U2 _ \  kh  y  (y2~y)  (Ti-T2) 

p l k i  P2k2  )  Hi-H2-cp2  ( T i -T 2 )  ( 121.1) 

1  a 


10* 


=  33,324. 


dynes 


cm 


assuming  a  =  .3961 
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The  sum  of  the  forces  is: 

+  17,020.  dynes 
cm2 

The  thermal  front  is  therefore  stable. 

However,  the  displacement  front  travels  faster  than  the  thermal 
front  and  is  therefore  isothermal.  The  displacement  front  stability 
criterion  is  simplified  from  Miller's  criterion  to  yield: 

(  .  pm  +  P2V2  \  108  „  _  (P2_P1)  g  a  >  0 

\  kl  k2  / 

For  Harmony  44  displacing  Harmony  87: 

+  ,25X,QQ56..  103  ( .3961 ) - (  .848-. 869)  981  ( .3961 ) 

100  100  ) 

=  -5098  d^nes 
cm7 

Therefore  the  displacement  front  is  unstable  at  experimental 
velocities. 


XII.  APPENDIX  F 


Raw  Data  and  Data  Work-Up 

This  section  contains  a  numerical  simulation  flow  chart,  an 
illustration  of  experimental  strategy,  examples  of  probability  plots  for 
experimental  and  numerical  results,  and  a  summary  of  raw  data.  Examples 
of  refractive  index  and  vicosity  calibration  curves  are  also  included. 
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Figure  29  Flowchart  of  Numerical  Program 
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t  =  0 


t  =  t 

+  At 

Itr  = 

Itr  +1 
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Figure  30  Experimental  Strategy 
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Figure  31  Probability  Plot  of  Dispersion  Profile,  R3 
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Figure  32  Probability  Plot  of  Numerical  Dispersion  Profile 
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CD 


Figure  33  Harmony  44~111,  Refractive  Index 
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Figure  34  varsol~Harmony  111,  Refr\  Index 
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Figure  35  Harmony  44-87,  Refractive  Index 
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Figure  36  Harmony  87~111,  Refractive  Index 
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Figure  37  Harmony  53~87  ,  Refractive  Index 
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Figure  38  Harmony  44~53  f  Refractive  Index 
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Figure  39  varsol“Harmony  44.  Refractive  Index 


EXPERIMENTAL 

LOGARITHMIC 

CRAGOE 


110 


Figure  40  Viscosity  of  f1ixtures,Harrnony  44~87 


TABLE  5 

Summary  of  Raw  Data,  Runs  1-36 
Run  1 

Harmony  111  displacing  Harmony  44 
TQ  =  26 . 7  C  T  1  =27 . 4  C 


T  i  me 

(  s ) 

Concentrat i on 
( %  Harmony  111) 

240. 

300. 

3960. 

4020. 

4080. 

4140. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5820. 

5880. 

5940. 

6000. 

.020 
.035 
.000 
.000 
.000 
.000 
.007 
.000 
.007 
.007 
.020 
.030 
.030 
.085 
.  135 
.206 
.268 
.341 
.422 
.575 
.699 
.752 
.834 
.915 
.932 
.957 
.970 
.985 
.995 
.990 
.985 
.995 

1  .000 

1  .0 

1  .00 

1  .00 

1  .0 

porosi ty= . 328  pore  vo1ume=1260  cm3 
80  -120  urn  sand  (20-30  mesh) 


' 
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Run  2 

Harmony  44  displacing  Harmony  111 
T0  =  27 . 5  T  1  =23 . 5 


T  i  me 

Concentration 

( s ) 

(%  Harmony  44) 

600. 

.576 

630. 

.757 

660. 

.548 

690. 

.576 

720. 

.565 

750. 

.630 

780. 

.585 

810. 

.570 

840. 

.550 

870. 

.527 

900. 

.522 

930. 

.613 

960. 

.545 

990. 

.580 

1020. 

.585 

1050. 

.576 

1080. 

.576 

1110. 

.522 

1140. 

.505 

1170. 

.522 

1200. 

.477 

1230. 

.465 

1260. 

.494 

1290. 

.494 

1320. 

.450 

1350. 

.438 

1380. 

.430 

1410. 

.387 

1440. 

.340 

1470. 

.325 

1500. 

.312 

1530. 

.315 

1560. 

.307 

1590. 

.290 

1620. 

.285 

1650. 

.285 

1680. 

.250 

1710. 

.240 

1740. 

.240 

1770. 

.250 

1800. 

.245 

1830. 

.240 

1860. 

.228 

1890. 

.240 

1910. 

.228 

1940. 

.240 

poros i ty 

=  .  328  pore  vo 

■ 

Run  3 

Harmony  111  displacing  Harmony  44 
Heated  injection,  T0  =  27  C  T 1  =92  C 


T  i  me 
( s ) 

Concentrat ion 
( %  Harmony  111) 

3600. 

4140. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5820. 

5880. 

.015 

.007 

.007 

.007 

.007 

.022 

.015 

.030 

.042 

.062 

.075 

.152 

.222 

.307 

.410 

.465 

.500 

.737 

.798 

.845 

.895 

.895 

.960 

.985 

.980 

.995 

.995 

1  .00 
.995 

1.00 

1.00 

Oil  flow  rate  =948  cm3/hr 

porosity  =.328  pore  volume  =1260  cm3 

80-120  urn  sand 


- 
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Run  4 

Harmony  44  displacing  Harmony  111 
T 0  =  2 5 . 3  C  T  1  =25 . 6  C 


T  i  me 

Concentrat ion 

( s ) 

(%  Harmony  44) 

180. 

1.00 

360. 

1.00 

540. 

.995 

720. 

1.00 

900. 

1.00 

1080. 

1.00 

1260. 

1  .00 

1440. 

.995 

1620. 

.990 

1800. 

.700 

1980. 

.748 

2169. 

.848 

2340. 

.740 

2520. 

.815 

2700. 

.715 

2880. 

.695 

3060. 

.730 

3240. 

.647 

3420. 

.630 

3600. 

.705 

3780. 

.630 

3960. 

.565 

4140. 

.565 

4320. 

.510 

4500. 

.515 

4680. 

.420 

4860. 

.420 

5040. 

.515 

5220. 

.460 

5400. 

.440 

5580. 

.365 

5760. 

.308 

5940. 

.308 

6120. 

.330 

6300. 

.302 

6480. 

.295 

6660. 

.273 

6840. 

.270 

7020. 

.245 

7200. 

.285 

porosi ty= 

.328  pore  volume 

80-120  urn 

i  sand  (20-30  mesh) 

Oil  flow 

rate  =937  cm3/hr 

1260  cm3 


' 
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Run  5 

Harmony  87  displacing  Harmony  44 
T0  =  27 . 6  C  T  1  =27 . 2  C 


T  i  me 

( s ) 

Concentration 
(%  Harmony  87) 

3600. 

4140. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5880. 

5940. 

.000 

.000 

.023 

.023 

.023 

.031 

.031 

.041 

.095 

.  135 

.  192 

.285 

.396 

.482 

.585 

.664 

.74  1 

.824 

.895 

.910 

.932 

.970 

.980 

.988 

.995 

.988 

.995 

.988 

1  .00 
.980 
.995 

porosity  =  .328  pore  volume  =  1260  cm3 
80-120  urn  sand  (20-30  mesh) 

Oil  flow  rate  =945  cm3/hr 

note : Concentr at i ons  were  normalized  to  1.0  and  0.0 


. 
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Run  6 

Harmony  111  displacing  Harmony  87 
T0=27 . 2  C  T 1 =26 . 7  C 


T  ime 

Concen 

( s ) 

(%  Har 

4560. 

.000 

4620. 

.000 

4680. 

.000 

4740. 

.052 

4800. 

.020 

4860. 

.020 

4920. 

.052 

4980. 

.052 

5040. 

.068 

5100. 

.115 

5160. 

.  130 

5220. 

.  145 

5280. 

.195 

5340. 

.305 

5400. 

.402 

5460. 

.450 

5520. 

.482 

5580. 

.465 

5640. 

.755 

5700. 

.755 

5760. 

.740 

5820. 

.917 

5880. 

.932 

5940. 

.968 

6000. 

.932 

6060. 

.987 

6120. 

.968 

6180. 

.987 

poros i ty= . 328  pore  volume=1260  cm3 

80-120  urn  sand 

Oil  flow  rate=851  cm3/hr 


- 

' 
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Run  7 

varsol  displacing  Harmony  111 
T0=27 . 0  C  T 1 =27 . 5  C 


T  ime 

Concentra 

( s ) 

(%  varsol 

180. 

.995 

360. 

.995 

7540. 

.995 

720. 

.945 

900. 

.596 

1080. 

.484 

1260. 

.436 

1440. 

.436 

1620. 

.447 

1800. 

.405 

1980. 

.328 

2160. 

.287 

2340. 

.255 

2520. 

.250 

2700. 

.243 

2880. 

.242 

3060. 

.232 

3240. 

.214 

3600. 

.213 

3780. 

.210 

3960. 

.202 

4140. 

.203 

4320. 

.204 

4500. 

.207 

4680. 

.  197 

5040. 

.  184 

5220. 

.  181 

5400. 

.  177 

5580. 

.  176 

5760. 

.  177 

5940. 

.  178 

6120. 

.  179 

6300. 

.  177 

6480. 

.  177 

6660. 

.177 

6840. 

.  178 

7020. 

.177 

7200. 

.  186 

poros i ty= . 328  pore  volume=1260.  cm3 
80-120  urn  sand 
Oil  flow  rate=970  cm3/hr 
Concentrations  were  normalized  to  1.0 


and  0 . 0 
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Run  8 

Harmony  111  displacing  varsol 
T0=28 . 2  C  T 1 =28  C 


T  i  me 

( s ) 

Concentration 
( %  Harmony  111) 

3660. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5820. 

5880. 

6000. 

.000 
.000 
.009 
.014 
.032 
.062 
.  1  17 
.210 
.362 
.486 
.604 
.778 
.771 
.910 
.952 
.957 
.976 
.985 
.999 
.997 
.997 
.997 
.997 
.998 

1.00 

1.00 

poros i ty= . 328  pore  volume=1260  cm3 
80-120  urn  sand 
Oil  flow  rate=963  cm3 /hr 
Concentrations  were  normalized  to  1.0  and  0.0 


■ 
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Run  9 

Harmony  87  displacing  Harmony  44 
Heated  injection,  T0  =  26  C  T 1 =  94  C 


T  ime 

( s ) 

Concent r at i on 
(%  Hamony  87) 

3680. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5940. 

.000 
.000 
.000 
.020 
.038 
.085 
.094 
.  175 
.265 
.352 
.442 
.575 
.622 
.782 
.830 
.855 
.925 
.940 
.972 
.990 
.995 
.990 

1  .00 

1.00 

.995 

1.00 

1  .00 

porosi ty= . 328  pore  volume=1260  cm3 

80-120  urn  sand 

Oil  flow  rate=980  cm3/hr 


Run  1  0 

Harmony  111  displacing  Harmony  87 
Heated  injection,  T0  =  26.5  C  T  1  =  9 6 . 0  C 


T  ime 

Conce 

(s) 

(%  Ha 

3660. 

.000 

4080. 

.000 

4140. 

.020 

4200. 

.000 

4260. 

.000 

4320. 

.020 

4380. 

.020 

4440. 

.020 

4500. 

.090 

4560. 

.  115 

4620. 

.  162 

4680. 

.210 

4740. 

.290 

4800. 

.385 

4860. 

.482 

4920. 

.500 

4980. 

.585 

5040. 

.688 

5100. 

.792 

5160. 

.810 

5220. 

.845 

5280. 

.898 

5340. 

.934 

5400. 

.950 

5460. 

.970 

5520. 

.987 

5580. 

1  .00 

5640. 

1  .00 

5700. 

1  .00 

6000. 

1.00 

poros i ty= . 328  pore  volume=1260  cm3 

80-120  urn  sand 

Oil  flow  rate=966  cm3/hr 
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Run  1  1 


Hamony  87 

displacing  H 

T0=26 . 1  C 

T 1 =26  C 

T  ime 

Concentrat i on 

( s ) 

(%  Hamony  87) 

14400. 

.000 

18000. 

.000 

18240. 

.000 

18480. 

.000 

18720. 

.012 

18960. 

.066 

19200. 

.076 

19440. 

.  103 

19680. 

.195 

19920. 

.285 

20160. 

.430 

20400. 

.528 

20640. 

.621 

20880. 

.753 

21120. 

.783 

21360. 

.845 

21600. 

.903 

21840. 

.924 

22080. 

.934 

22320. 

.980 

22560. 

1  .00 

22800. 

1  .00 

23040. 

1.00 

23760. 

1  .00 

poros i ty= 

. 328  pore  vo 

80-120  um 

sand 

low  flow 

rate 

44 


vol ume  = 1260 


cm3 


' 


Run  1 2 

Hamony  111  displacing  Harmony  87 
Chilled  injection,  T  0  =  80  C  T  1  =  1 9 . 5  C 


T  ime 

Concentrat i on 

is) 

( %  Harmony  111) 

3360. 

.000 

4200. 

.000 

4260. 

.000 

4320. 

.000 

4380. 

.017 

4440. 

.033 

4500. 

.069 

4560. 

.  158 

4620. 

.230 

4680. 

.318 

4740. 

.389 

4800. 

.578 

4860. 

.578 

4920. 

.653 

4980. 

.789 

5040. 

.808 

5100. 

.925 

5160. 

.944 

5220. 

.967 

5340. 

.985 

5400. 

.985 

5460. 

1  .00 

5520. 

1  .00 

5580. 

1  .00 

6000. 

1  .00 

poros i ty= . 328  pore  volume=1260  cm3 
80-120  urn  sand 

Concentrations  were  normalized  to  1.0  and  0.0 


Run  1  3 

Hamony  44  displacing  Harmony  111 
Heated  injection,  T  0  =  26.0  C  T 1 =  9 0  C 


Time  Concentration 


(  s  )  (  % 

1260. 

1440. 

1620. 

1800. 

1980. 

2160. 

2340. 

2520. 

2700. 

2880. 

3060. 

3240. 

3420. 

3600. 

3780. 

3960. 

4140. 

4320. 

4500. 

4680. 

4860. 

5040. 

5220. 

5400. 

5580. 

5760. 

5940. 

6120. 

6300. 

6480. 

6660. 

6840. 

7020. 

7200. 


Harmony  44) 

005 

005 

150 

170 

315 

310 

160 

125 

252 

315 

420 

451 

516 

505 

545 

522 

545 

590 

607 

607 

652 
658 

675 
660 

653 
653 

676 
698 
704 
698 
704 
698 
738 
732 


porosi ty= . 328  pore  volume=1260  cm3 

80-120  urn  sand 

Oil  flow  rate=979  cm3/hr 


. 

' 

Run  1  4 

Harmony  87  displacing  Harmony  44 
T0  =  26 . 3  C  T 1 =26 . 7  C 


T  ime 

Conce 

( s ) 

(%  Ha 

3660. 

.000 

4080. 

.000 

4140. 

.000 

4200. 

.000 

4260. 

.005 

4320. 

.012 

4380. 

.046 

4440. 

.100 

4500. 

.  193 

4560. 

.282 

4620. 

.360 

4680. 

.493 

4740. 

.640 

4800. 

.692 

4860. 

.815 

4920. 

.870 

4980. 

.894 

5040. 

.915 

5100. 

.963 

5160. 

.970 

5220. 

.988 

5280. 

.988 

5340. 

.988 

5400. 

1.00 

5460. 

1.00 

5520. 

.996 

5580. 

1  .00 

5640. 

1  .00 

6000. 

1.00 

porosi ty= . 328  pore  volume=1260  cm3 

80-120  urn  sand 

Oil  flow  rate=975  cm3 /hr 


. 
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Run  1 5 

Harmony  111  displacing  Harmony  87 
T0  =  28 . 0  C  T  1  =28 . 0  C 


T  i  me 

(s) 

Concentration 
( %  Harmony  111) 

3660. 

3840. 

3900. 

3960. 

4020. 

4080. 

4140. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4720. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

6000. 

.000 
.000 
.000 
.000 
.007 
.038 
.022 
.038 
.053 
.  100 
.117 
.  133 
.230 
.297 
.297 
.424 
.475 
.634 
.634 
.791 
.808 
.808 
.841 
.881 
.936 
.936 
.963 
.982 

1  .00 
.982 

1  .00 

1  .00 

1.00 

poros i ty= . 328  pore  volume=1260  cm3 

80-120  urn  sand 

Oil  flow  rate=964  cm3/hr 

concentrations  were  normalized  to  1.0  and  0.0 


' 
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Run  1 6 

Harmony  111  displacing  Harmony  44 
T0=26 . 9  C  T 1 =27 . 8  C 


T  ime 

( s ) 

Concentrat ion 
( %  Harmony  111) 

3660. 

4140. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5940. 

.000 

.000 

.000 

.000 

.007 

.007 

.048 

.068 

.110 

.195 

.250 

.370 

.494 

.602 

.685 

.778 

.845 

.880 

.900 

.942 

.942 

.950 

.962 

.962 

.985 

.985 

.985 

.995 

poros i ty= . 328  pore  volume=1260  cm3 

80-120  urn  sand 

Oil  flow  rate=963  cm3/hr 


' 

Run  1  7 

Hamony  111  displacing  Harmony  87 
T0  =  27  C  T  1=27.2  C 


T  i  me 

Concent 

( s ) 

(%  Harm 

3660. 

.000 

4320. 

.000 

4380. 

.000 

4440. 

.000 

4500. 

.017 

4560. 

.017 

4620. 

.051 

4680. 

.  180 

4740. 

.218 

4800. 

.273 

4860. 

.419 

4920. 

.419 

4980. 

.557 

5040. 

.597 

5100. 

.735 

5160. 

.752 

5220. 

.784 

5280. 

.816 

5340. 

.834 

5400. 

.834 

5460. 

.957 

5520. 

.957 

5580. 

.957 

5640. 

.957 

5700. 

.977 

5820. 

1.000 

5880. 

1  .000 

6000. 

1  .000 

poros i ty= , 328  pore  volume=1260  cm3 
80-120  urn  sand 
oil  flow  rate=927 
Concentration  was  normalized  to  1.0 


cm3/hr 
norma  1 i zed 


and  0.0 


Run  1 8 

Harmony  87  displacing  Harmony  44 
T0  =  27 . 0  C  T  1=26.8  C 


T  ime 
( s ) 

Concentration 
(%  Harmony  87) 

3660. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5820. 

5880. 

5940. 

6000. 

6060. 

6120. 

6180. 

6240. 

6300. 

.000 

.000 

.000 

.000 

.010 

.032 

.090 

.175 

.302 

.446 

.597 

.680 

.835 

.934 

.934 

.970 

.993 

1  .00 
.993 

1  .00 

1  .00 

1.00 

1.00 

poros i ty= . 368  pore  volume=1411  cm3 
20-30  urn  (80-120  mesh)  sand 
oil  flow  rate=908  cme/hr 


. 
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Run  1 9 

Harmony  111  displacing  Harmony  87 
T  0  =  24.9  C  T  1  =25 . 3  C 


T  i  me 

Conce 

( s ) 

(%  Ha 

4860. 

.035 

5100. 

.000 

5400. 

.020 

5460. 

.085 

5520. 

.085 

5580. 

.130 

5640. 

.  146 

5700. 

.257 

5760. 

.322 

5820. 

.450 

5880. 

.567 

5940. 

.620 

6000. 

.755 

6060. 

.810 

6120. 

.810 

6180. 

.933 

6240. 

.950 

6300. 

.950 

6360. 

.950 

6420. 

.970 

poros i ty= . 368  pore  vo1ume=1411  cm3 
20-30  urn  sand 


RUN  20 


Harmony  44  displacing  1 

o 

CM 

CM 

II 

O 

1— 

11=24.4  C 

T  i  me 

Concent rat i on 

( s ) 

(%  Hamony  44) 

00180. 

.005 

00360. 

.005 

01980. 

.005 

02160. 

.005 

02340. 

.150 

02520. 

.  105 

02700. 

.212 

2880. 

.160 

03060. 

.  165 

03240. 

.  195 

03420. 

.212 

03600. 

.222 

03780. 

.222 

03960. 

.295 

04140. 

.360 

04320. 

.370 

04500. 

.420 

04680. 

.425 

04860. 

.435 

05040. 

.470 

05220. 

.480 

05400. 

.535 

05580. 

.545 

05760. 

.590 

05940. 

.590 

06120. 

.630 

06300. 

.635 

06480. 

.635 

06660. 

.640 

06840. 

.675 

7020. 

.698 

7200. 

.705 

poros i ty= 

. 368  pore  vo 

20-30  urn 

sand 

Oil  flow 

rate=900  cm3/ 

1  cm3 


' 

Run  21 

Harmony  87  displacing  Harmony  44 
Heated  injection,  T0  =  25.5  C  T  1=91  .5  C 


T  ime 

( s ) 

Concent r at i on 
(%  Harmony  87) 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5880. 

.000 

.000 

.022 

.085 

.156 

.246 

.407 

.556 

.832 

.903 

.960 

.994 

1.00 

1.000 

1  .00 

1  .00 

poros i ty= . 368  pore  volume=1411  cm3 

20-30  urn  sand 

oil  flow  rate=972  cm3/hr 


. 
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Run  22 

Harmony  111  displacing  Harmony  87 
Heated  injection,  T 0  =  28 -  1  C  T 1=92.4  C 


T  i  me 

(s) 

Concentration 
( %  Harmony  111) 

4260. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5820. 

5880. 

5940. 

6000. 

6600. 

.000 
.000 
.000 
.024 
.000 
.054 
.  134 
.  134 
.  167 
.284 
.365 
.479 
.619 
.691 
.778 
.816 
.835 
.871 
.928 
.962 
.962 
.962 
.979 

1.00 

1.00 

1.00 

poros i  ty=  .  368  pore  volume=141.1  cm3 

20-30  urn  sand 

oil  flow  rate=950  cm3/hr 

Concentrations  were  normalized  to  1.0  and  0.0 
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Run  23 

Harmony  44  displacing  Harmony  111 
Heated  injection,  T0=28.7  C  T 1=89.0  C 


Concentration 
(%  Harmony  44) 


T  i  me 
( s ) 

0180. 

01440. 

1620. 

1800. 

1980. 

2160. 

2340. 

2520. 

2700. 

2880. 

3060. 

3240. 

3420. 

3600. 

3780. 

3960. 

4140. 

4320. 

4500. 

4680. 

4860. 

5040. 

5220. 

5400. 

5580. 

5760. 

5940. 

6120. 

6300. 

6480. 

6660. 

6840. 

7020. 

7200. 


.000 
.000 
.000 
.047 
.091 
.136 
.  145 
.  179 
.205 
.205 
.224 
.273 
.310 
.354 
.364 
.354 
.389 
.414 
.467 
.480 
.530 
.563 
.575 
.574 
.581 
.586 
.591 
.636 
.662 
.662 
.690 
.702 
.729 
.747 


poros i ty= . 368  pore  volume=1411  cm3 

20-30  urn  sand 

Oil  flow  rate=925  cm3/hr 

Concentrat ions  were  normalized  to  1.0  and  0.0 


Run  24 

Harmony  87  displacing  Harmony  44 
Chilled  injection,  T0=80.3  C  T  1  =26 . 9  C 


T  ime 
( s ) 

Concentrat i on 
(%  Harmony  87) 

4860. 

5100. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5820. 

5880. 

5940. 

6000. 

6060. 

6120. 

6180. 

6240. 

6300. 

6360. 

6600. 

.046 
.000 
.046 
.046 
.067 
.074 
.  142 
.181 
.315 
.450 
.600 
.756 
.837 
.920 
.936 
.936 
.982 

1  .00 

1  .00 

1  .00 

1.00 

poros i ty= . 368  pore  volume=1411  cm3 

20-30  urn  sand 

oil  flow  rate=855  cm3/hr 

Concentrations  were  normalized  to  1.0  and  0.0 


Run  25 

Harmony  87  displacing  Harmony  44 
T0  =  30 . 0  C  T  1  =30 . 0  C 


T  i  me 

( s ) 

Concentrat ion 
(%  Harmony  87) 

21600. 

23040. 

23280. 

23540. 

23760. 

24000. 

24240. 

24480. 

24720. 

24960. 

25200. 

25440. 

25680. 

25920. 

26180. 

26420. 

26660. 

26900. 

27140. 

27380. 

27620. 

28340. 

.000 
.000 
.000 
.011 
.000 
.011 
.045 
.  108 
.  190 
.31  1 
.405 
.494 
.602 
.670 
.761 
.847 
.892 
.926 
.943 

1  .00 

1  .00 

1  .00 

poros i ty= . 368  pore  volume=1411  cm3 

20-30  urn  sand 

Low  flow  rate  (210  cm3/hr) 

Concentrations  were  normalized  to  1.0  and  0.0 


. 
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Run  26 

Harmony  87  displacing  Harmony  44 
T0  =  29  C  T  1  =29 . 3  C 


T  i  me 


Concentrat ion 


(s) 

( %  Har 

65700. 

.000 

72900. 

.000 

73800. 

.000 

74700. 

.000 

75600. 

.000 

76500. 

.011 

77400. 

.022 

78300. 

.072 

79200. 

.  172 

80100. 

.250 

81000. 

.327 

81900. 

.544 

82800. 

.588 

83700. 

.672 

84600. 

.758 

85500. 

.844 

86400. 

.917 

87300. 

.950 

88200. 

1  .  00 

89100. 

1  .00 

90000. 

1  .00 

poros i ty= . 368  pore  volume=1411  cm3 
20-30  urn  sand 

Field  displacement  rate  (61  cm3/hr) 
Concentrations  were  normalized  to  1 


and  0 . 0 


' 


137 


Run  27 

water  wet  sand  pack 

Harmony  87  displacing  Harmony  44 

T0=29 . C  T 1 =29 . 2  C 


Concentration 
( %  Harmony  87 ) 


T  ime 

( s ) 

3660. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5460. 


.000 
.000 
.000 
.028 
.055 
.075 
.  170 
.274 
.432 
.545 
.692 
.787 
.871 
.921 
.936 
.952 
.952 
1  .00 
.995 
.995 
1  .00 


poros i ty= . 356  pore  vol.=1354  cm3 
volume  of  residual  water  in  pack=174  cm3 
oil  flow  rate  =  975  cm3/hr 
Oil  concentrations  are  normalized  to  0.0  andl.O 


' 


Run  28 

Harmony  111  displacing  Hamony  87 

water  wet  pack  20-30  um  sand  (80-120  mesh) 

T0  =  27 . 2  C  T 1  =27 . 5  C 


Concentration 
( %  Harmony  111) 


T  ime 

( s ) 

3660. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5940. 


.000 
.000 
.000 
.043 
.043 
.  116 
.  189 
.226 
.371 
.429 
.560 
.597 
.719 
.777 
.899 
.859 
.940 
.959 
1.00 
1  .00 
1.00 
1  .00 


poros i ty= . 356  pore  vol.=1354  cm3 
volume  of  residual  water=75  cm3 
oil  flow  rate=975  cm3/hr 

Concentrations  are  normalized  to  1.0  and  0.0 


. 
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Run  29 

water  wet  sand  pack 

Harmony  44  displacing  Harmony  111 

T0=27 . 5  C  T 1 =27 . 8  C 


T  ime 
( s ) 

Concent rat i on 
(%  Harmony  44) 

360. 

540. 

900. 

1260. 

1620. 

1800. 

1980. 

2160. 

2340. 

2520. 

2700. 

2880. 

3060. 

3400. 

3780. 

4320. 

4680. 

5040. 

6840. 

7020. 

.000 

.000 

.000 

.000 

.000 

.005 

.039 

.157 

.222 

.264 

.292 

.365 

.519 

.481 

.562 

.632 

.686 

.697 

.768 

.773 

poros i ty= . 356  pore  vol ume= 1 354cm3 

The  volume  of  residual  water  could  not  be  calculated 
oil  flow  rate=975  cm3/hr 

Concentrations  were  normalized  to  1.0  at  endpoint. 


Run  30 

Heated  displacement,  T0  =  26.  C  T 1 =  92 . 0  C 
Harmony  44  displacing  varsol 


Time  Concentration 

(s)  (%  Harmony  44) 


3660. 

3720. 

4080. 

4140. 

4200. 

4260. 

4320. 

4380. 

4440. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

6000. 


.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.018 
.026 
.041 
.  1  16 
.223 
.331 
.420 
.570 
.770 
.883 
.931 
.958 
.970 
.990 
1  .00 
1  .00 
1  .00 
1.00 


porosi ty= . 356  pore  volume=1354.  cm3 
20-30  urn  sand  (80-120  mesh) 

Some  aqueous  phase  was  present  in  the  pack. 
Calculated  volume=27.  cm3 
The  concentr at i ons  of  the  end  points 
were  normalized  to  1.0  and  0.0 


Run  3  1 

Harmony  87  displacing  Harmony  44 
T0  =  25 . 2  C  T  1  =25 . 8  C 


T  i  me 

(s) 

Concent rat i on 
(%  Harmony  87) 

3600. 

3900. 

4200. 

4500. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5700. 

6000. 

.000 

.000 

.000 

.000 

.000 

.020 

.020 

.080 

.145 

.145 

.370 

.535 

.690 

.810 

.880 

.920 

1.00 

1.00 

1  .00 

1  .00 

poros i ty= . 356  pore  volume=1354  cm3 
oil  flow  rate=970.  cm3/hr 
20-30  urn  sand 


' 
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Run  32 

Harmony  111  displacing  Harmony  87 
T0  =  27 . 0  C  T  1  =27 . 5  C 


Time  Concentration 

( s )  ( %  Harmony  111) 


3660. 

3900. 

4200. 

4500. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5520. 

5580. 

5640. 

5700. 

5760. 

5820. 

5880. 

5940. 

6060. 


.025 
.000 
.000 
.000 
.000 
-.056 
.000 
.018 
.037 
.  194 
.233 
.342 
.420 
.582 
.624 
.726 
.810 
.828 
.873 
.873 
.936 
.958 
.978 
1.00 
1.00 


poros i ty= . 356  pore  volume=1354  cm3 
20-30  urn  sand  (80-120  mesh) 
oil  flow  rate  =940  cm3/hr 

Concentrations  were  normalized  to  0.0  and  1.0 


at  endpoints. 


. 
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Run  33 

Harmony  111  displacing  Harmony  44 

Heated  displacement  T0=27  C  T1=86.4  C 


T  i  me 

Conce 

(s) 

(%  Ha 

4740. 

.000 

4800. 

.000 

4860. 

.016 

4920. 

.036 

4980. 

.095 

5040. 

.213 

5100. 

.323 

5160. 

.462 

5220. 

.546 

5280. 

.753 

5340. 

.829 

5400. 

.930 

5460. 

.943 

5520. 

.963 

5580. 

.950 

5640. 

.985 

5700. 

.985 

5760. 

1.00 

5820. 

1.00 

5940. 

.995 

6000. 

1.00 

Oil  flow  rate  =975  cm3 /hr 
Porosity  =.356  pore  volume=1354  cm3 
20-30  urn  (80-120  mesh)  sand 
concentrations  were  normalized  to  1.0  and  0.0 


at  end  points. 


. 

Run  34 

Harmony  53  displacing  Harmony  44 
Heated  injection,  T0  =  25  C  T  1=92.5  C 


Time  Concentration 

(s)  (%  Harmony  53) 


4200. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5700. 

6000. 


.000 
.000 
.000 
.025 
.025 
.069 
.  167 
.243 
.389 
.531 
.645 
.772 
.861 
.890 
.861 
.984 
.984 
.984 
1  .00 
1  .00 


oil  flow  rate  =1020  cm3/hr 
porosity  =.356  pore  volume=1354 
20-30  urn  sand 


cm3 


■ 


Run  35 

Harmony  87  displacing  Harmony  53 
Heated  displacement,  T0=25  C  T1 


T  i  me 

Conce 

(s) 

(%  Ha 

3660. 

.000 

3900. 

.000 

4200. 

.000 

4500. 

.000 

4800. 

.000 

4860. 

.010 

4920. 

.076 

4980. 

.  123 

5040. 

.244 

5100. 

.408 

5160. 

.532 

5220. 

.720 

5280. 

.846 

5340. 

.871 

5400. 

.888 

5460. 

.990 

5520. 

.990 

5580. 

1  .00 

5700. 

1.00 

6000. 

1  .00 

89.2  C 


poros i ty= . 356  pore  volume=1354  cm3 
20-30  urn  (80-120  mesh)  sand 


' 


Run  36 

Harmony  44  displacing  varsol 
Heated  injection,  T0  =  25  C  T 1 =  9 1 .5  C 


T  ime 

( s ) 

Concentrat ion 
(%  Harmony  44) 

3660. 

3900. 

4200. 

4500. 

4560. 

4620. 

4680. 

4740. 

4800. 

4860. 

4920. 

4980. 

5040. 

5100. 

5160. 

5220. 

5280. 

5340. 

5400. 

5460. 

5640. 

6000. 

.000 
.000 
.000 
.005 
.015 
.020 
.031 
.066 
.  153 
.306 
.463 
.648 
.814 
.890 
.949 
.969 
.995 
.995 

1  .00 

1.00 

1  .00 

1.00 

pore  volume=1354  cm3  poros i ty= . 356 
20-30  urn  sand  (80-120  mesh) 


